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Abstract 

We study weakly stable semilinear hyperbolic boundary value problems with highly oscil- 
latory data. Here weak stability means that exponentially growing modes are absent, but the 
so-called uniform Lopatinskii condition fails at some boundary frequency j3 in the hyperbolic 
^ region. As a consequence of this degeneracy there is an amplification phenomenon: outgoing 

■^ , waves of amplitude O(e^) and wavelength e give rise to reflected waves of amplitude 0(e), so the 

r^ ' overall solution has amplitude 0(e). Moreover, the reflecting waves emanate from a radiating 

wave that propagates in the boundary along a characteristic of the Lopatinskii determinant. 

An approximate solution that displays the qualitative behavior just described is constructed 
by solving suitable profile equations that exhibit a loss of derivatives, so we solve the profile 
equations by a Nash-Moser iteration. The exact solution is constructed by solving an associated 
singular problem involving singular derivatives of the form dx' + /3-^, x' being the tangential 
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variables with respect to the boundary. Tame estimates for the linearization of that problem 
are proved using a first-order calculus of singular pseudodifferential operators constructed in 
T^lj- ' the companion article |CGW12] . These estimates exhibit a loss of one singular derivative and 

C^ ■ force us to construct the exact solution by a separate Nash-Moser iteration. The same estimates 

#y\ ' are used in the error analysis, which shows that the exact and approximate solutions are close 

f^ . in L°° on a fixed time interval independent of the (small) wavelength e. The approach using 

OQ ' singular systems allows us to avoid constructing high order expansions and making small divisor 

assumptions. 
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1 Introduction and main results 

In this paper we study weakly stable semilinear hyperbolic boundary value problems with oscil- 
latory data. The problems are weakly stable in the sense that exponentially growing modes are 
absent, but the uniform Lopatinskii condition fails at a boundary frequency /3 in the hyperbolic 
region ?^. As a consequence of this degeneracy in the boundary conditions, there is an amplifica- 
tion phenomenon: boundary data of wavelength e and amplitude 0{e'^) in the problem ()l.ip below 
gives rise to a response of amplitude 0{e). In the meantime, resonance may occur between distinct 
oscillations. In the situation studied below, a resonant quadratic interaction between two incoming 
waves of amplitude 0{e) may produce an outgoing wave of amplitude O(e^). When reflected and 
amplified on the boundary, this oscillation gives rise to incoming waves of amplitude 0(e). Hence 
the 0(e) amplitude regime appears as the natural weakly nonlinear regime. 

Let us now introduce some notation. On M_,_ = {x = {x',Xd) = {t,y,Xd) = {t,x") : Xd > 0}, 
consider the N x N semilinear hyperbolic boundary problem for v = V£{x), where e > Qj: 

(a) Lo{d)v + fo{v) = 0, 

(1.1) ^i>)cl>{v)=e^G(x',^) onxd = 0, 

(c) V = and G = m t < 0, 



^See Definition 11.31 and Assumption 1 1.51 below for precise statements. 
We usually suppress the subscript e. 



where LQ{d) = dt + ^j=iBj dj, the matrix B^ is invertible, and both fo{v) and (/)(f) vanish at 
V = 0. The function G{x',6o) is assumed to be periodic in ^O; and the frequency /? G M^ \ {0} is 
taken to be a boundary frequency at which the so-called uniform Lopatinskii condition fails. A 
consequence of this failure is that the choice of the factor e^ in (jl.ip (b) corresponds to the weakly 
nonlinear regime for this problem. The leading profile is nonlinear ly coupled to the next order 
profile in the nonlinear system (jl.35p . (|1.36p derived below. We also refer to Appendix [B] for a 
detailed specific example which illustrates the nonlinear feature of the leading profile equation. 

Before proceeding we write the problem in an equivalent form that is better adapted to the 
boundary. After multiplying (jl.ip (a) by [Bd)~^ we obtain 

L{d)v + f{v) = 0, 

(1.2) ^{v)=e^G(x',^-^^ oiixd = 0, 



V = and G = in t < 0, 



where we have set 



d-l 

L{d) = dd + Y^ Aj dj with Aj := B^^Bj for j = 0, . . . , d - 1. 

3=0 

Setting V = eu and writing f{v) = D{v)v, (j){v) = ip{v)v, we get the problem for u = Us{x) 

(a) L{d)u + D{eu)u = 0, 

(1.3) (6) V(en)u = eG fx', ^^ J on Xd = 0, 

(c) u = in t < 0. 

For the problem (II. 3p we pose the two basic questions of rigorous nonlinear geometric optics: 

(1) Does an exact solution u^ of (jl.3p exist for e G (0,1] on a fixed time interval [0,To] 
independent of e? 

(2) Suppose the answer to the first question is yes. If we let n"^^ denote an approximate 
solution on [0, Tq] constructed by the methods of nonlinear geometric optics (that is, solving eikonal 
equations for phases and suitable transport equations for profiles) , how well does U£-^^ approximate 
Ue for e small? For example, is it true thaCl 

(1.4) lim|ue-Mf^|Loo^O? 

The amplification phenomenon was studied in a formal way for several different quasilinear 
problems in )AM871 [MA881 IMR83] . In |MR83] amplification was studied in connection with Mach 
stem formation in reacting shock fronts, while |AM87j explored a connection to the formation 
of instabilities in compressible vortex sheets. Both papers derived equations for profiles using 
an ansatz that exhibited amplification; however, neither of the two questions posed above were 



Let us observe that by the amphfication phenomenon, we expect the solution v to (II. ip to have ampUtude 0{e), 
so the solution u to (|1.3|l should have amplitude 0(1). Hence the limit (|1.4p deals with the difference between two 
0(1) quantities. 



addressed. The first rigorous amplification results were proved in [CGlOj for linear problems. That 
article provided positive answers to the above questions (question (1) is trivial for linear problems) 
by making use of approximate solutions of high order, and showed in particular that the limit (|1.4p 
holds. 

In this paper we give positive answers to the above questions for the semilinear system p.3p . 
As is typical in nonlinear geometric optics problems involving several phases, difficulties with small 
divisors rule out the construction of high order approximate solutionci- Instead of constructing the 
exact solution u^ as a small perturbation of a high-order approximate solution, we construct u^ in 
the form 

Ue{x) = Ue{x,6o)\^^^p^, 

where U£{x,9o) is an exact solution of the singular system (jl.lSp . The singular system is solved 
using symmetrization and diagonalization arguments from |Wil02] , modified and supplemented with 
methods from |Cou04j for deriving linear estimates for weakly stable hyperbolic boundaryproblems. 
In deriving the basic estimate (12. 4p for the singular linear problem, a loss of derivativeqj forces us 
to use a new tool, namely, a substantial refinement, given in the companion paper |CGW12] . of 
the calculus of singular pseudodifferential operators constructed in |Wil02j . In the new version of 
the calculus, residual operators have better smoothing properties than previously realized and can 
therefore be considered as remainders in our problem. The loss of derivatives in the linear estimate 
presents a serious difficulty in the application to our semilinear problem. Picard iteration appears 
to be out of the question, so in section 15.21 we use a Nash-Moser iteration scheme adapted to the 
scale of spaces ()1.19p to construct the solution Us{x,6q) to the semilinear singular problem. 

If the problem (jl.3p satisfied the uniform Lopatinskii condition, then because of the factor e 
in the boundary data eG, the equations for the leading profile, V*^ in (ll.lSp . would be linear; and 
in fact V'^ would vanish. The weakly nonlinear regime would correspond to a source term G (and 
not eG) in ()1.3p . see |Wil96l fWilO O] . Under our weak stability assumption, it turns out that V*^ is 
nonlinearly coupled to the second-order profile V^ in the profile equations (jl.35p . (|1.36p . To solve 
these equations we first isolate a "key subsystem" (jl.42p that decouples from the full system. The 
basic L^ estimate for the linearization of the key subsystem still exhibits a loss of one derivative, 
and we are again forced to use Nash-Moser iteration in order to solve this subsystem. Once the key 
subsystem is solved, the solution of the full profile system p.35p . (ll.36p follows easily. It appears 
in our analysis that the leading order amplitude equation shares the weak well-posedness of the 
original nonlinear problem but we have not checked whether the loss of derivative for the amplitude 
equation is optimal (we conjecture it is so). 

The error analysis used to answer question (2) above is based on the estimate for the singular 
system (jl.lSp . see Proposition 12.21 below, and is discussed in more detail in section [T31 

This paper can be read independently of [CGW12] : for the reader's convenience, we have gath- 
ered all the necessary material on the singular calculus in Appendix [Aj Before giving a fuller 
discussion, we first provide some definitions, notation, and a precise statement of assumptions. 

1.1 Assumptions 

We make the following hyperbolicity assumption on the system (jl.ip : 



*Such difficulties are sometimes avoided by assuming that small divisors do not occur, see e.g. [JMR93) . but we 
do not want to make this assumption. 

^In fact, the basic L^ estimate for the singular system H1.18I) exhibits loss of a single "singular derivative" d^i H , 

which is optimal according to the analysis in [CGlOj . 



Assumption 1.1. There exist an integer g > 1, some real functions Xi, . . . , Xq that are analytic on 
M.'^\ {0} and homogeneous of degree 1, and there exist some positive integers vi, . . . ,Vq such that: 

d q 

VC = (6, . . . ,ed) G IK' \ {0} , det [r/ + ^e, Bj\ = H (^ + ^'^(^))'' ■ 

j=i fc=i 

Moreover the eigenvalues Ai(^), . . . , Aq(^) are semi- simple (their algebraic multiplicity equals their 
geometric multiplicity) and satisfy Ai(^) < • • • < Ag(^) for all ^ ^W^\ {0}. 

For simplicity, we restrict our analysis to noncharacteristic boundaries and therefore make the 
following: 

Assumption 1.2. The matrix B^ is invertihle and the matrix B := ip{0) has maximal rank, its 
rank p being equal to the number of positive eigenvalues of Ba (counted with their multiplicity). 
Moreover, the integer p satisfies 1 < p < N — 1. 



In the normal modes analysis for (jl.3p . one first performs a Laplace transform in the time 
variable t and a Fourier transform in the tangential space variables y. We let r — 27 G C and 
T] G M denote the dual variables of t and y. We introduce the symbol 



d-l 

-1 



A{0■.= -iB^'\{T-i^)I + Y,l^,BJ\ , C:=(r-i7,^)GCx 

For future use, we also define the following sets of frequencies: 

H := {(r - i7,7?) G C X M'^-i \ (0, 0) : 7 > o} , S := {c G H : r^ + 7^ + |r?|2 = l} , 

Ho := { (r, 7?) G M X R^'^ \ (0, 0)} = H n {7 = 0} , So : 



sn: 



Two key objects in our analysis are the hyperbolic region and the glancing set that are defined 
as follows: 

Definition 1.3. • The hyperbolic region % is the set of all {T,rf) G Hq such that the matrix 
A{t, Tj) is diagonalizable with purely imaginary eigenvalues. 

• Let G denote the set of all (r, ^) G M x M'^ such that ^ 7^ and there exists an integer 
k G {1, . . . , g} satisfying: 

// 7r(G) denotes the projection of G on the d first coordinates (in other words ■7t{t, S,) = 
(r, ^1, . . . ,S,d-i) for all {t,S,)), the glancing set Q is Q := vr(G) C Ho- 

We recall the following result that is due to Kreiss |Kre70j in the strictly hyperbolic case (when all 
integers Vj in Assumption 11.11 equal 1) and to Metivier |MetOO| in our more general framework: 

Proposition 1.4 ( |Kre701 IMetOO| ). Let Assumptions \1.1\ and \l.'^ be satisfied. Then for all Q G 
H \ Ho, the matrix A{C,) has no purely imaginary eigenvalue and its stable subspace E*(C) has 
dimension p. Furthermore, W defines an analytic vector bundle over H \ Hq that can be extended 
as a continuous vector bundle over H. 



For all {T,r]) G Hq, we let W{T,r]) denote the continuous extension of E** to the point [t^t]). The 
analysis in [MetOO j shows that away from the glancing set ^ C Sq, IE*((^) depends analytically on 
Q, and the hyperbolic region % does not contain any glancing point. 

To treat the case when the boundary operator in (jl.3p (b) is independent of n, meaning ip{eu) = 
'(/'(O) =: B, we make the following weak stability assumption on the problem {L{d),B). 

Assumption 1.5. • For allC (^^\ ^o, KerS n E'*(C) = {0}. 

• The set To := {C G So : Keri? n E*(C) 7^ {0}} is nonempty and included in the hyperbolic 
region %. 

• For all C, G Tq, there exists a neighborhood V of ( inT,, a real valued C°° function a defined 
on V, a basis Ei{(^), . . . , Ep{C,) ofW{C,) that is of class C°° with respect to ^ £ V, and a matrix 
-P(C) £ GLp(C) that is of class C°° with respect to (^ £ V, such that 

VCGV, B {E.iO ... EpiC))=P{Odmg{j + ia{C),l,...,l). 

For comparison and later reference we recall the following definition. 

Definition 1.6 ( |Kre70| ). As before let p be the number of positive eigenvalues of B^. The problem 
{L{d),B) is said to be uniformly stable or to satisfy the uniform Lopatinskii condition if 

B : E"(C) -^CP 

is an isomorphism for all ^ £ S. 

Remark 1.7. Observe that if {L{d),B) satisfies the uniform Lopatinskii condition, continuity 
implies that this condition still holds for {L{d),B+ip), where ip is any sufficiently small perturbation 
of B. Hence the uniform Lopatinskii condition is a convenient framework for nonlinear perturbation. 
The analogous statement may not be true when (L{d),B) is only weakly stable. Remarkably, 
weak stability persists under perturbation in the so-called WR class exhibited in [BGRSZ02], and 
Assumption [L5] above is a convenient equivalent definition of the WR class (see |CG10^ Appendix 
B]). In order to handle general nonlinear boundary conditions as in (jl.3p we shall strengthen 
Assumption 11.51 in Assumption 11.111 below. 

Boundary and interior phases. We consider a planar real phase (po defined on the boundary: 

(1.5) 4>oit,y) ■.= Tt + ri-y, (r,7?)GHo. 

As follows from earlier works, see e.g. [MA88J . oscillations on the boundary associated with the 
phase <j)o give rise to oscillations in the interior associated with some planar phases (pm- These phases 
are characteristic for the hyperbolic operator -Lo(^) ^'iid their trace on the boundary {x^ = 0} equals 
(j)Q. For now we make the following: 

Assumption 1.8. The phase (J)q defined by (jl.5p satisfies {T,r]) G To- In particular {T,rj) G Ti. 

Thanks to Assumption ll.81 we know that the matrix A{t, rj) is diagonalizable with purely imaginary 
eigenvalues. These eigenvalues are denoted iuii, . . . i^^Mi where the w^'s are real and pairwise 
distinct. The w^'s are the roots (and all the roots are real) of the dispersion relation: 

d-l 

det zl + ^'n_Bj+ujBd =0. 



To each root w^ there corresponds a unique integer /c^ S {1, . . . ,(/} such that r + Afe,^(7],a;„) = 0. 
We can then define the following reajj phases and their associated group velocities: 

(1.6) Vm = l,...,M, (t)m{x) :=Mt:y)+klmXd, v^ := VAfe^(ry,u;„) . 

Let us observe that each group velocity v^ is either incoming or outgoing with respect to the space 
domain M*^: the last coordinate of v^ is nonzero. This property holds because {t,!]) does not 
belong to the glancing set Q. We can therefore adopt the following classification: 

Definition 1.9. The phase (p^ is incoming if the group velocity v^ is incoming (that is, when 
'^id^kmi'Hi^m) > ^J) '^^'^ ^^ ^-5 outgoing if the group velocity v^ is outgoing (d^^Xk^{r],ui^) < 0). 

In all that follows, we let I denote the set of indices m £ {1, . . . , M} such that (pm is an incoming 
phase, and O denote the set of indices m G {1, . . . , M} such that 4'm is an outgoing phase. If p > 1, 
then I is nonempty, while if p < A^ — 1, O is nonempty (see Lemma ll.lOl below). We will use the 
notation: 

d d-l 

j=l k=0 

P-={Z,V), x' = {t,y), 4>o{x') = ^ ■ x' . 

For each phase (pm, dp'm denotes the differential of the function 0^, with respect to its argument 
X = {t,y,X(i). It follows from Assumption 11.11 that the eigenspace of A{(3) associated with the 
eigenvalue iw^ coincides with the kernel of Loldcpm) and has dimension z^^^^^. The following well- 
known lemma, whose proof is recalled in |CG10| . gives a useful decomposition of W in the hyperbolic 
region. 

Lemma 1.10. The stable subspace W{(3) admits the decomposition: 

(1.7) E^(/3) = e™exkerLo(d</.„), 

and each vector space in the decomposition (|1.7p admits a basis of real vectors. 

To formulate our last assumption we observe first that for every point C, £% there is a neigh- 
borhood V of ^ in S and a C°° conjugator Qq{C) defined on V such that 

Awi(C)/ni \ 

(1.8) Qo(C)^(C)Qo '(0 = ••. =:-Di(C), 

V iu^AOlnj) 

where the ujj are real when 7 = and there is a constant c > such that either 

Re (iujj) < —cy or Re (iuj) > 07 for all C £ 1^- 

In view of Lemma ll. 101 we can choose the first p columns of Qq (C) to be a basis of E*((^), and we 
write 

Qo\0 = [Q^n{OQout{0]■ 



^If (Hi v) does not belong to the hyperbolic region H, some of the phases (pm niay be complex, see e.g. [Wil96l 
IWilOOl [LesOTI [MarlO] . Moreover, glancing phases introduce a new scale ^/e as well as boundary layers. 

7 



Choose J' so that the first J' blocks of — Bi he in the first p columns, and the remaining blocks in 
the remaining N — p columns. Thus, Re (iwj) < —cy for 1 < j < J'. 
Observing that the linearization of the boundary condition in (|1.3|) is 

u I — > ip{eu)u + [dvip{£u) ii] eu , 

we define the operator 

(1.9) B{vi,V2)u := il){vi)u + [dv^{vi)u]v2 , 
which appears in Assumption 11.11] below. For later use we also define 

(1.10) :P(i^i, 1^2) ii := D{vi)u + [d^D{vi) ii] V2 , 
as well as 

(1.11) B{vi) := B{vi,vi) , V{vi):=V{vi,vi). 

We now state the weak stability assumption that we make when considering the general case 
of nonlinear boundary conditions in (jl.Sp . 



Assumption 1.11. • There exists a neighborhood O o/(0,0) G M^^ such that for all (^1,^2) G 
O and allQ G H \ Hq, ker 6(^1,^2) n E^(C) = {0}. For each (^1,^2) G O the set T{vi,V2) := 
{C G So : ker B{vi,V2) n ]E*(C) ^ {0}} is nonempty and included in the hyperbolic region %. 
Moreover, if we set T := ^{v-i,v2)eo'^{'^ii'^'2), then T dl-L (closure in SqJ. 

• For every (" G T there exists a neighborhood V of C in T, and a C°° function cr{vi,V2,C) on 
O xV such that for all (t^i, i'2, () G O x V we have Ker B{vi,V2) n E*(C) / {0} if and only if 
C G So CLiid a{vi,V2,C) = 0. 

Moreover, there exist matrices Pi(wi,W2,C) £ GLp(C), i = 1,2, of class C°° on O x V 
such that V(fi,U2,C) G O x V 

(1.12) Pi{vi,V2,OB{vi,V2)Qin{OP2{vi,V2,C) = diag {-f + ia{vi,V2,0A,---A) ■ 

For nonlinear boundary conditions, the phase (po in (|1.5p is assumed to satisfy (r, ?]) G T(0, 0), 
or in other words the intersection ker i? n E* (r , ry) is not reduced to {0} (the set Tq in Assumption 
11.51 is a short notation for T(0, 0)). The phases (j)m are still defined by (|1.6p and thus only depend 
on L{d) and B and not on the nonlinear perturbations /q and ip{£u) — V'(O) added in (II. 3p . 

Remark 1.12. 1) The properties stated in Assumption 1 1 . 1 ij are just a convenient description of the 
requirements for belonging to the WR class of |BGRSZ02] . Like the uniform Lopatinskii condition. 
Assumption 1 1.11] can in practice be verified by hand via a "constant-coefficient" computation. More 
precisely, for (ui,U2) near (0,0) G M^^ and C S S, one can define (see, e.g., [BGS071 chapter 4]) a 
Lopatinskii determinant A(z;i, V2, C) that is C°° in (tii, ^2)) analytic in (^ = (r — 27, 77) on T,\Q, and 
satisfies 

A(ui,';;2,C) = Oif andonly if Kere(ui,i;2)nE''(C) / {0}. 

In particular, A(t;i,f2, •) is real-analytic on Ti. 



Following |BGRSZ02] , see also |BGS071 chapter 8] , we claim that Assumption 11.111 holds pro- 
vided 

(1.13) 0/{CG S: A(0,0,C) =0} cn and A(0,0,C) = 0^a^A(0,0,C) ^0, 

and it thus only involves a weak stability property for the linearized problem at (1^1,^2) = (0,0). 
Indeed, the implicit function theorem then implies that for (ui,i'2) near zero and {T,r]) near (", the 
set {(r, 77) G Eq : A{vi,V2,t,7]) = 0} is a real-analytic hypersurface in Ti. On the other hand, an 
application of the implicit function theorem to A{vi, V2, z,r]), for (z, 77) G E shows that the real 
dimension of the manifold {{z, ry) G S : A(fi, ■U2, z, ry) = 0} must be the same, that is, d—2. The two 
zero sets must then coincide; there are no zeros in S \ Sq. The function a and the neighborhoods O 
and V arise in a factorization of A given by the Weierstrass Preparation Theorem. The construction 
of the conjugating matrices Pi, i = 1,2 follows from a construction in |ST88l p. 268-270]. 

Instead of assuming (11.130 , we have stated Assumption 11.111 in a form that is more directly 
applicable to the proof of Proposition 12.21 and to the error analysis of Theorem 14.11 

2) To prove the basic estimate for the linearized singular system. Proposition 12. 2^ and to 
construct the exact solution Us to the singular system (J1.18p below, it is enough to require that the 
analogue of Assumption 1 1 . 1 1 1 holds when B{vi,V2) is replaced hy B{vi) := B{vi,vi). However, for 
the error analysis of section [J] in the case of nonlinear boundary conditions, we need Assumption 
11.111 as stated. 

The next Lemma, proved in [CGlOj . gives a useful decomposition of C^ and introduces projec- 
tors needed later for formulating and solving the profile equations. 

Lemma 1.13. The space C^ admits the decomposition: 

(1.14) C^ = e^iikerLo(#,„) 

and each vector space in ()1.14p admits a basis of real vectors. If we let Pi,- ■ ■ , Pm denote the 
projectors associated with the decomposition (|1.14p . then there holds Im B^ Lo{d4)m) = kerP^ for 
allm = l,...,M. 

1.2 Main results 

For each m G {1, . . . , M} we let 

^m,ky rC — i, . . . , I'km 

denote a basis of ker LQ{d(l)m) consisting of real vectors. In section[5]we shall construct a "corrected" 
approximate solution u^ of (II. 3p of the form 

(1.15) utix) = V' U^\ +eV' U^\ +e'u'Jx,^,^Y 
where (p := (01, . . . , (J)m) denotes the collection of all phases, 

^ ^ mex fc=i ^ ^ 



(1.16) 



m=l k=l 



and the crm,k{x,6m) and Tm,k{x,6m) are scalar C^ functions periodic in 6m with mean which 
describe the propagation of oscihations with phase (pm and group velocity Vm- Here TZ denotes the 
operator 

7^V° = -R [L(5^)V° + L»(0)V°] 

for R defined as in ()1.32p . The last corrector e^Wp(x,6'o,Cd) iii (|1-15|) is a trigonometric polynomial 
constructed in the error analysis of section [H 

The following theorem, our main result, is an immediate corollary of the more precise Theorem 
KT\ Here we let Qt ■= {{x, 6*0) = {t, y, Xd, Oq) G M'^+i x T^ : x^ > 0, t < T}, and bQr ■= {{t, y, Oq) S 
M'^ X T^ : t < T}. The spaces E^ are defined in IHTidh below. 

Theorem 1.14. We make Assumptions HOI IJ.^l \1.5[ and \1.8\ when the boundary condition in 
(|1.3|) is linear (ip{£u) = ip{0)); in the general case we substitute Assumption \l.lI\ for Assumption 
[731 Fix T >0, set Mq := 3 d + 5, and let 



a := [{d + l)/2] + M0 + 3 andd = 2a- [{d + l)/2]. 

Consider the semilinear boundary problem (jl.Sp . where G{t,y,6Q) £ H'^[hQ,T)- There exists £q> Q 
such that if {G)ija+2n,Q_^\ is small enough, there exists a unique function Ueix^Oo) G E°'~^{Qt) 
satisfying the singular system (ll.lSp on ily such that 

Ue[x) := Ue X, 



is an exact solution of (II. 3p on (— oo,T] x M^ for < e < eg- -^'^ addition there exists a profile 
V^{x,9) as in ()1.16p . whose components cjrn,k Us in H'^~^{Q,t), such that the approximate solution 
defined by 



uf P := V 



satisfies 

lim \ue - <pp|loo =0 on (-00, T] x m1. 

Observe that although the boundary data in the problem (jl.3p is of size 0(e), the approximate 
solution Ue^^ is of size 0(1), exhibiting an amplification due to the weak stability at frequency 
/3. The main information provided by Theorem 11.141 is that this amplification does not rule out 
existence of a smooth solution on a fixed time interval, that is it does not trigger a violent instability, 
at least in this weakly nonlinear regime. As far as we know, the derivation of the leading order 
amplitude equation (jl.42p is also new in the general framework that we consider. We hope that the 
analysis developed in this article will be useful in justifying quasilinear amplification phenomena 
such as the Mach stems or kink modes formation of [AM871 IMA88] . 

Remark 1.15. a) In order to avoid some technicalities we have stated our main result for a problem 
(jl.3p where all data is in t < 0. This result easily implies a similar result in which outgoing waves 
defined in t < of amplitude 0(e) and wavelength e give rise to reflected waves of amplitude 
0(1). In either formulation, analysis of the profile equations (see Remark I1.27P shows that the 
waves of amplitude 0(1) emanate from a radiating wave that propagates in the boundary along a 
characteristic of the Lopatinskii determinant. 
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b) We have decided to fix T > at the start and choose data small enough so that a solution 
to the nonlinear problem exists up to time T. One can also (as discussed in Remark 13 .Tp fix the 
data in the problem {G in (|1.3p ) at the start, and then choose T small enough so that a solution 
to the nonlinear problem exists up to time T. 

In the remainder of this introduction, we discuss the construction of exact solutions, the con- 
struction of the approximate solution V'', and the error analysis. Complete proofs are given in 
Sections [21 [3l H and O 

1.3 Exact solutions and singular systems 

The theory of weakly stable hyperbolic initial boundary value problems fails to provide a solution 
of the system (jl.Sp that exists on a fixed time interval independent of e|II- In order to obtain such an 
exact solution to the system (jl.Sp we adopt the strategy of studying an associated singular problem 
first used by |JMR95] for an initial value problem in free space. We look for a solution of the form 

(1-17) Uejx) = [/e(x,6'o)|^^_ ^„(^/) , 

where Us{x,9q) is periodic in ^o and satisfies the singular system derived by substituting ()1.17p 
into the problem (|1.3p . Recalling that L{d) = dd + X]i=o ^j^j ^^ obtain: 

(1-18) ddU, + A (d,, + ^] Ue + D{eUe) U, = 

i'{eUe)Ue\^,=o = eG{x',eo), 
Us = Omt< 0. 

The special difficulties presented by such singular problems when there is a boundary are described 
in detail in the introductions to |Wil02] and [CGWll] . In particular we mention: (a) symmetry 
assumptions on the matrices Bj appearing in the problem p.ip equivalent to ()1.3p are generally of 
no help in obtaining an L^ estimate for (|1.18p (boundary conditions satisfying Assumption II . 51 can 
not be maximally dissipative, see |CG10j ): (b) one cannot control L°° norms just by estimating 
tangential derivatives d?,^ ^C/e because (ll.lSj) is not a hyperbolic problem in the Xd directioio; 
moreover, even if one has estimates of tangential derivatives uniform with respect to e, because of 
the factors 1/e in (ll.lSp one cannot just use the equation to control ddUs and thereby control L°° 
norms. 

In |Wil02j a class of singular pseudodifferential operators, acting on functions U{x' , Oq) periodic 
in 6q and having the form 

poU{x'M = —^Y.I e^"'-^'+*'^'''>(^l^(^',^o),C' + — ,7) t/(^',A;)dC', 7 > 1, 



^This would be true even for problems {L{d), B) that are uniformly stable in the sense of Definition [LB] 
®For initial value problems in free space, one can control L°° norms just by estimating enough derivatives tangent 
to time slices t — c. 
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was introduced to deal with these difficulties. Observe that the differential operator A appear- 
ing in (jl.lSp can be expressed in this form. Kreiss-type symmetrizers rs{Dx'^eo) in the singular 
calculus were constructed in [Wil02| for (quasilinear systems similar to) (|1.18|) under the assump- 
tion that (L(9),'0(O)) is uniformly stable in the sense of Definition 1 1.6i With these one can prove 
L^{xd, H^{x' , 9q)) estimates uniform in e for (|1.18|) . even when e G is replaced by G in the boundary 
condition. To progress further and control L°° norms, the boundary frequency (3 is restricted to 
lie the complement of the glancing set. With this extra assumption the singular calculus was used 

in [WHO 2] to block-diagonalize the operator A (eU£,dx' H — ) microlocally near the /3 direction 

and thereby prove estimates uniform with respect to e in the spaces 

(1.19) E' ■= C{xd, H%x', 9o)) n L\xa, H'+\x', Oo)). 

These spaces are Banach algebras and are contained in L°° for s > — ^. For large enough s, as 
determined by the requirements of the calculus, existence of solutions to (ll.lSp in E^ on a time 
interval [0, T] independent of e E (0, Eq] follows by Picard iteration in the uniformly stable case. 

The singular calculus of |Wil02j was used again in [CGWl lJ to rigorously justify leading order 
geometric optics expansions for the quasilinear analogue of (II. 3p in the uniformly stable case (with 
(3 £ Ti and the forcing term G in place of eG in the boundary condition). Under the assumptions 
made in the present paper, in particular assuming weak stability as in Assumption 11.51 or Assump- 
tion [TTTTl we face the additional difficulty that the basic L^ estimate for the problem {L{d),B) 
exhibits a loss of derivatives. A consequence of this is that the singular calculus of [Wil02j is no 
longer adequate for estimating solutions of (jl.lSp . The main reason is that remainders in the 
calculus of [Wil02] are just bounded operators on L^, while for energy estimates with a loss of 
derivative remainders should be smoothing operators. We therefore need to use an improved ver- 
sion of the calculus constructed in |CGW12] in which residual operators are shown to have better 
smoothing properties than previously thought. With the improved calculus we are able in section 
12.31 to estimate solutions of (jl.lSp in E^ spaces (|1.19p , but of course there is a loss of one singular 
derivative in the estimates. This loss forces us in section 15.21 to use Nash-Moser iteration on the 
scale of E^ spaces to obtain an exact solution of the singular system (jl.lSp on a fixed time interval 
independent of e. Observe that one singular derivative costs a factor 1/e and this is another reason 
why the scaling eG in (jl.lSp is crucial. 

Remark 1.16. The main idea for proving the estimate for the linearized singular problem. Propo- 
sition [221 is to adapt the techniques of |Cou04] to the singular pseudodifferential framework. There 
is however one major obstacle along the way. While the error term in the composition of two zero 
order operators (or in the composition of an operator of order of order —1 (on the left) with an 
operator of order 1, a (—1, 1) composition) is smoothing of order one in the sense of (lA.Sp . the same 
is unfortunately not true of the error term in (1, —1) compositions (there are counter-examples for 
that). The properties of the (1, —1) error terms that arise in our proof are described in ()2.32p : see 
also Proposition IA.9I To deal with these errors we use further microlocal cutoffs x*^ in the extended 
calculuqj and partition the solution as in (12. 6p . In particular we are led to estimate terms like those 
involving norms of V^j/A^ ^2 in ™ Proposition 12.61 and V^^/A^ U^ in Proposition! 



These are cutoffs in (^',m)-space which depend on (e,7); they are described in section lA. 41 
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1.4 Derivation of the leading profile equations 

We now derive the profile equations for tlie semilinear problem (|1.3|) . We work with profiles 
V^ {x, 6) periodic in 6 = (6*1, . . . , 6m), where 9j is a placeholder for (j)j/e. Looking for an approximate 
solution of (|1.3p of the form u"" = (V^ + eV^ + e'^V'^)\g=^/^, where <j) = (01, . . . , (pAi), we get interior 
equations 

(a) £(9e)V° = 0, 

(1.20) (6) C{de)V^ + L(a)V° + D(0)V° = 0, 

(c) C{de)V^ + L{d)V^ + D{0)V^ + (a^L>(0)V°)V° = 0, 

by plugging u"" into (|1.3p (a) and setting the coefficients of, respectively, e~^, e°, and s equal to 
zero. The operator C^de) is defined by 

M 

(1.21) £{de):=J2L{d(t>j)dey 

With B := tp{0), the boundary equations, obtained by plugging u"" into (ll.3j) (b) and setting the 
coefficients of e^ and e equal to zero, are 

i3V°(x',0,eo,...,^o) = 0, 
^ ' ^ i3Vi + (a,V(O)VO)V° = G(x',0o), 

where ^o is a placeholder for (po/^- We will see that as a consequence of the weak stability at 
frequency /3, the problem for the leading profile V*^ is nonlinear and nonlocal. (See Appendix [B] for 
a concrete example.) Thus, the scaling in (11. 2p is the weakly nonlinear scaling when the uniform 
Lopatinskii condition fails at a hyperbolic frequency /3. To analyze these equations we proceed to 
define appropriate function spaces and a pair of auxiliary operators E and R. 
Functions V{x,6) G L'^(R^ x T*^) have Fourier series 

(1.23) Vix,9)= Y, Va{x)e'^-'. 

Since only quadratic interactions appear in (ll.20p and we anticipate that V*^ will have the form in 
(fTTG]) . for A; = 1,2 we let 

Z ' = {a G Z : at most k components of a are nonzero}, 
and we consider the subspace i7*;^(M+^^ x T^^) C i^^(l+^^ x T^^) defined by 



(1.24) H'-^''{R'^^^ xT^^)= h{x,e)£H'{R'^^'^ xT^'^):V{x,e)= J^ Va{x)e'' 



Thus, multiplication defines a continous map 

(1.25) ii'^'i(l++^ X T^^) X /7^'i(M++^ X T^) ^ F"'2^1'J+^ X T^) 

for s > (d + 1 + 2)/2. 
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Definition 1.17. Setting (j) := ((/>!, . . . ,(J)m), we say a G Z '^ is a characteristic mode and write 
a gC if det L(d{a ■ (j))) = 0. Otherwise we call a a noncharacteristic mode. We decompose 

C = U^^^Cm, where Cm := {a G Z ''^ : a ■ (j) = nacf)^ for some n^ G Z}. 

Observe that for a G Cm, the integer n^ is necessarily equal to Ylik=i^k- Since (j)i and 4)j are 
linearly independent for i ^ j^ any a G Z^"^'^ \ belongs to at most one of the sets C^ and n-Q, 7^ 
if a / 0. 

Elements a G Cm with two nonzero components correspond to resonances. Resonances are 
generated in products like crp^k{x, ^)o'r,k'{x, ^), which arise from the quadratic term in ()1.20p (c). 
whenever there exists a relation of the form 

T^m4'm = T^p4'p + n^(t>r^ where m G {1, . . . , M} \ {p, r} and n^, n^, n^ G Z. 

We then refer to (</>m) (t>p^ (pr) as a triple of resonant phases. This relation implies, for example, 
that (pp oscillations interact with (pr oscillations to produce (pm oscillations. 

Definition 1.18. We define the continuous projecto^ E : H''^'^(W^\^^ x T^^) -^ H''^(^^^^ x T^^), 
s>Q,hy 

M 

(1.26) E = Eq+^ Em, where EqV -.= ¥0 and EraV := ^ PmVa(x)e*""'^'", 

m=l aeCm\0 

for Pm as in Lemma \1.13X 

For C{de) as in (fL2n) we have that for V° G iJ'^'2(R'J_+^ x T*0, 

(1.27) ^V° = V° if and only if V° G ^^^^(mJ^^ x T*^) and C{de)V^ = 0, 
and ()1.27p in turn is equivalent to the property that V'' has an expansion of the form 

(1.28) V° = v{x) +^Y^ ^"^''^ (^' ^'"^ ""^'^ ' 

m=l A;=l 

for some real- valued functions c^m^k- Moreover, since for any ttt., 

d-l 

(1.29) L{dcPm) =ui^I + Y,l^3 ^0 = E (^- - ^'fc)^'^' 

we have for V G iJ'^^^^^J+i ^ y^^): 

(1.30) ^/:(5e)V = C{de)EV = 0. 

We also need to introduce a partial inverse R for Lido). We begin by defining 



°Thc continuity of E is shown in [CGWll] , Remark 2.5. 
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which in view of (jl.29p satisfies 

(1.31) L{d4)m)Rm = RmLidcpm) = I - Pm- 



rrf+1 



The operator R is defined formally at first on functions V{x, 9) = Ylaez'^'-'^ Va{x)e^°' of //*' (M^ x 
T^) by 

(1.32) RV:= Yl R{a)Vaix)e'''-^ 



a& 



where 



(1.33) R{a) := < 



' lh^"i if a G Cm \ {0}, M 

if a = 0, and C{ia) := i >^ L{d(j)j)aj = iL{d{a ■ (j))). 

£{ia)-^ if a ^ C, i=^ 



Remark 1.19. The operator R is well-defined on functions V S H''''^{R_^_ x T*^) whose spectrum 
contains only finitely many noncharacteristic modes, and then RV lies in the same space. Oth- 
erwise, there can be a problem with small divisors; the possibility of there being infinitely many 
noncharacteristic modes a for which detL(d(a • </>)) is close to zero can prevent convergence of 
(fL32]l in iJ*;2(Mj^^ x T^^ for any t. 

It follows readily from (fL3T]l that for J" € H''^{R'^^^ x T*^), s > 0, 

(1.34) C{de)RT = RC{de)T = (/ - E)T. 

Such J^ have no noncharacteristic modes. Along with (ll.30p equation (I1.34P implies 

Proposition 1.20. Suppose T € iJ"'HK+^^ x T^^), s > 0. Then the equation C{dg)V = F has a 
solution V G i7^'i(l+^^ X T^'^) if and only if EF = 0. 

By applying the operators E and R to the equations (|1.20p and using (|1.27p , (|1.30p , and (jl.34p , 
we obtain: 

(a) EV° = V°, 

(6) E{L{d)V^ + D(0)V°) = 0, 
(c) SV° = on Xd = 0, = (00 
{d) V° = in t < 0. 



(1.35) 



; • • • ) "\^J^ 



and 



(1.36) 



(a) (/ - E)V^ + R{L{d)V^ + D{0)V^) = 0, 

(b) E {L{d)V^ + D{0)V^ + {d^D{0) V°)V°) = 0, 

(c) BV^ + (a,V'(0)V°)V° = Gonxd = 0,e = {Bo, . . . , ^o), 

(d) V^ = in t < 0. 
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Remark 1.21. (a) Since EV^ = V° the function L{d)V° + D{0)V^ in (fL36l) fa') has no noncharac- 
teristic modes so the action of R on this function is weh-defined. 

(b) It is easy to check that functions V°, V^ belonging to H'^'^'^{R^ x T*^), s > ^, and sat- 
isfying (fL35]) and (fOBD fa) also satisfy ([00]) (a), (b) and ([1:22]) . Equation [L36]^b) and Proposition 
11.201 suggest that we might obtain a solution of ()1.20|) (c) by taking 

(/ - E)V^ = -R (L{d)V^ + D{0)V^ + (a^Z)(0)V°)V° 

There are two problems with this. First, the quadratic term {dyD{0)V^)V^ generally has infinitely 
many noncharacteristic modes, so one should expect a problem with small divisors. Second, the 
statement (jl.34p and Proposition 11.201 are both not true when J^ G H^''^(R_^_ x T ), even if J^ 
has finitely many noncharacteristic modeq^j. These difficulties affect the error analysis and are 
discussed further in section 11.51 

To determine the equations satisfied by the individual profiles v{x), i7m,fc(a^, ^m) in the ex- 
pansion (11.28P of V'^, we first refine the decomposition of the projector E in (11.26p . For each 
m G {1,... ,M} we let 

denote a basis of real vectors for the left eigenspace of the real matrix 

d-l 

(1.37) a(/3) = rAo + ^7?.^, 



■J ^ 



associated to the eigenvalue — w^, chosen to satisfy 

\ 1, ii m = m' and k = k', 

^m,k ' '^m' ,k' \ , 

lU, otherwise. 

For V G C^ set 

Pm,kV ■■= iim,k ■ v)rm,k (no complex conjugation here). 

We can now write 



E — Eq + 2_^ 2^ Em,k, 



m=l k=l 

where -Em,fc '■= Pm,kEra- When the multiplicity /c = 1 we write E^a instead of E^a,! and do similarly 

for (.m,k, rm,k and so on. 

The following lemma, which is a slight variation on a well-known result [Lax57j , is included for 
the sake of completeness: 

Lemma 1.22. Suppose EV^ = V^ and that V^ has the expansion (11.28p . Then 

Em,k{L{d)V'^) = {X^^am,k)rm,k 



^^This is because of the fact that for any k £ Z\ {0}, there can be many a G {Cm \ 0) n Z'^''^ such that Ua = k. 
See the proof of Proposition 11.281 
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where X^^ is the characteristic vector field associated to ^mLj; 

d-l 
j=0 

Proof. For ^' £ Ti near /3, let —ujm{C') be the eigenvalues iA{^') ()1.37p and Pm{£,') the correspond- 
ing projectors (these objects depend smoothly on ^' near /3 thanks to the analysis of (MetOOj ) . 
Differentiate the equation 

d-l \ 

with respect to ^j, evaluate at /?, apply P^ on the left to obtain 

from which the lemma readily follows. D 

By Assumption 11.51 we know that the vector space keri? n W{/3) is one-dimensional; moreover, 
it admits a real basis because B has real coefficients and ]E*(/3) has a real basis. This vector space 
is therefore spanned by a vector e G M^ \ {0} that we can decompose in a unique way by using 
Lemma ll.lOl 

(1.38) ker B D E'(/3) = Span {e} , e = ^ e„ , P^ e^ = e,^ . 

Each vector e^ in (|1.38p has real components. We also know that the vector space BW{I3) is 
(p — l)-dimensional. We can therefore write it as the kernel of a real linear form: 

(1.39) BE'{^) = {X eCP, b-X = 0}, 

for a suitable vector b £W\ {0}. 

Any function V{x,9) G H^''^{M.^ x T ) can be decomposed 

V = V + Vine + Vout + Vnonch = V + V*, 

where the terms correspond respectively to the parts of the Fourier series (I1.23P with a = 0, a 
incoming, a outgoing, and a noncharacteristico. 

Proposition 1.23. Suppose V° G H''''^{M.'^^^ x T^^), s>l, is a solution of (fTB^ . Then 

V° = 0, VL = 0, Vl„,, = 0, andsoV'^ = VL = EVL, 

V {x' , 0,6q, . . . , 9q) = a{x' , 9q) e for some unknown periodic function a with mean 0. 



^^The vector field X^^ is a constant multiple of the vector field dt + Vm ■ ^x" computed by Lax for the Cauchy 
problem, where Vm is the group velocity defined in Definition 11.91 
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Here we say a is incoming if a G Cm \ for an index m such that 0™, is an incoming phase. 



17 



Proof. Since EV^ = V'', we have V^^^j^ = 0. Applying Eq to the problem (jl.35p . we find that the 
mean value V satisfies the weakly stable boundary problem 

L(a)V° + Z)(0)V° = 
BV° = on Xd = 
V° = in t < 0. 

By the well-posedness result of |Cou05j we have V*^ = 0. 

Lemma [1.221 implies that outgoing profiles am,k, m G O, in the expansion (jl.28p of V" satisfy 
problems of the form 

"km 

fe'=l 

crm,k = in t < 0, 

where X,^^ is an outgoing vector field. Thus, am,k = for all /c = 1, . . . , Vk^- 

Part (b) follows immediately from the boundary condition in (I1.35P and (jl.38p . D 

Since V° = V^^^, we obtain from (jl.36p (a) 

(/ - E)V' = {I- E)VL = -R{Lid)V'> + D{0)V'), 

so 

V'=y}+ Vi, + VL G H''\ where EV^, = V^- 

Next decompose the boundary condition ()1.36p (c): 

(1.40) BEVL = G*- [(a,^(0)V°)V°)]* - BVLt - B{I - E)vL 

= G*- [(a„V'(0)VO)V°)]* - BV^^t + BRiLid)V^ + D(0)V°) . 

Remark 1.24. (a) If V^uf\xa=o,e =eo were known, one could write down a transport equation for 
a{x',6o) which is determined by the solvability condition for (|1.4Up implied by (|1.39p : 

(1.41) b-{G* - [(a^V(0)V°)V°)]* - BVl^t + BR{L{d)V^ + D(0)V°)) = 0. 

However, the presence of the term E {^{dvD{0)V^)V^) in (]1.36p (b) implies that two incoming modes 
in Vf„^ (which is still unknown) can resonate to produce an outgoing mode that will affect V^^^f. 
Thus, we do not know Vluf\xa=o,ej=eo^ aiid we see that the nonlinear boundary equation (I1.41J) is 
coupled to the nonlinear interior equation ()1.36p . 

(b) If the phases are such that an outgoing mode can never be produced by a product of two 
incoming modes, then V^^^ can be determined from ()1.36p to be 0, and one can proceed as in [CGIO] 
to solve for a without having to use Nash-Moser iteration. 

The key subsystem to focus on now is (recalling V'' = EV^ = V^^^ and writing with obvious 
notation E = Eq + Einc + Eout)- 



(1.42) 



(a)i?,„,(L(a)VL + ^(0)VD = 0, 

(6) Eout {L{d)Vlt + DiO)VL, + (a,D(o)vL)vL) = o, 

(c) 6 • (G* - [(a,^(0)V°)V°)]* - BVl, + BR{L{d)Vl, + ^(0)VL)) = 0, 

(d)VL(x',O,0o,...,^o) = a(x',eo)e, 
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where V^^^ and V^^^ both vanish in t < 0. 

A formula for V^^^ in terms of a{x' , 60) can be determined by solving transport equations using 
(|1.42p (a). and that formula can be plugged into (|1.4'2|) (b) to get V^^j in terms of a. Thus, the 
subsystem ()1.42p can be expressed as a very complicated nonlinear, nonlocal equation for the 
single unknown a. This is done in Appendix [B] for a strictly hyperbolic example with only one 
resonance. However, that is not the way we solve (jl.42p : instead we solve the subsystem in its 
given form by iteration. Picard iteration does not work; there is a loss of derivatives from one 
iterate to the next (because of R), so we use a Nash-Moser scheme. An essential point is to 
take advantage of the smoothing property of the interaction integrals that pick out resonances in 
Eout{{dvD{0)Vf^^)V^^jLJi, that property allows us to get tame estimates in section [3l 

An important tool in solving the subsystem (jl.42p is the following result from [CGIO], which 
will allow us to write the boundary equation (jl.42p (c) as a transport equation for a{x',0Q). 

Proposition 1.25 ([CGIO], Proposition 3.5). Let the vectors b and e^ be as in (|1.39p . (jl.38p . and 
let (t(C) be the function appearing in Assuniption \1.5[ Then there exists a nonzero real number k 
such that 

R^Prn = for all m £ {1, ... , M} 

b ■ B 2^ RmAoCm = KdrO'{T, 7]) and drcriz, rj) = 1 , 

b- B'^ RmAjCm = Kdrj^a{T,ri), j = l,...,d-l. 



m£X 



and thus 



/ d-i \ 

b- B^ RmL{d)em = k dra{z,v)dt + ^5^,cr(r,r/)9x, =: Xlop- 

mgX \ j=l I 

Taking note of the factor -. — in the definition ()1.33p of R, we obtain the immediate corollary: 
Corollary 1.26. The boundary term b ■ B RL{d)V^^^ in (jl.42p may be written 

b-BRL{d)Vtc = XLopA, 
where A{x' , 9q) is the unique function with mean in 6q such that dg^A = a. 

Remark 1.27. Proposition 11.25) shows that propagation in the boundary, which is described by 
a{x' ,0o), is governed by the (x-projection of the) Hamiltonian vector field associated to the Lopatin- 
skii determinant. Since V''(rc', 0, ^O; • • • j^o) = CLix',Oo), this shows that waves of amplitude 0(1) 
emanate from the radiating boundary wave defined by a. 

After (I1.42P is solved V*^ is known, so V^, V^^^, and (/ — E)V}^^ can now be determined by 
returning to the full system (jl.36p . The trace of EV}^^ is not yet determined; one should make 
a choice of EVlji^\x^=o,ej=eo such that (jl.40p holds, and then solve for EVl^^ using (ll.36p (b). A 
precise description of the regularity of V" and V^ is given in Theorem 15.111 The last piece of the 
corrected approximate solution, e^U'i in (|1.15p is discussed next. 



'interaction integrals are similar to convolution integrals. 
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1.5 Error analysis 

Given a periodic function f{x, 9), where 9 = {9i, . . . , 9m), let us denote 

fix, 9)\g_,^g^^^^-j := fix, 9q + Wi^d, . . . , 6^0 + iilMid)] 

so we have 



Taking the profiles V*^, V^ constructed in Theorem 15. Ill if we define 

Ulix,9^) := iV''ix,9) + eV\x,9)) le^(,,,f^), 

we find that U^ satisfies the singular system 

(a) heiU^) := ddU^ + A (d^, + ^"j U^^ + DieU^,) Ul 
^^•^^^ ih) i)ieUl)Ul = eG(x' ,9q) + 0(e^) on x^ = 






(c) ^/° = in t < 0, 

where the error terms refer to norms in E^ and H^ spaces whose orders are made precise in section 
m For example, (jl.43p follows directly from the profile equations (|1.20p (a).(b). together with the 
identity 

(1.44) 

Le (/(x,0)|,^(,„,^)) = ^ iCide)fix,9)) |,^(,^^^) + iLid)fix,9)) |,_(,^^^) + p(e/)/)|,^(,„,^). 

Since our basic estimate for the linearized singular system exhibits a loss of one singular deriva- 
tive (basically, we lose a 1/e factor), the accuracy in (jl.43p (a) is not good enough to conclude that 
\Us — ^elL°o(x.eo) 1^ small (the error terms are only 0(e)). To improve the accuracy we construct 
an additional corrector Upix,9o,^d) and replace U^ by 



(1.45) Ueix, 9o) := (V°(x, 9) + eV\x, 9)) le^(eo,£^) + ^^^pi^, Oo, —)■ 



In constructing lA^ we deal with the first (small divisor) problem described in Remark ll.21l fbl by 
approximating V and V^ by trigonometric polynomials Vp and Vi to within an accuracy 5 > 
in appropriate Sobolev norms, and seek U'i in the form of a trigonometric polynomial^. To deal 
with the second (solvability) problem, we use the following Proposition, which allows us to use the 
profile equation (jl.36p (b) as a solvability condition, in spite of the failure of Proposition 1 1 . 201 when 

Mdeo,d^d) ■= Lid^o)deo + d^^. 



^^Trigonometric polynomial approximations were already used to deal with small divisor problems in the error 



analysis of [JMR95] . 
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Proposition 1.28. Suppose F{x,6) G H^''^{W_^_ x T ) has a Fourier series which is a finite sum 
and that EF = 0. Then there exists a solution of the equation 

(1-46) Co{de„d^,Mx,eo,^d) = F{x,e)\e^^eo,i,) 

in the form, of a trigonometric polynomial in (^O)'^d) of the form 

(1.47) U{x,eo,^d)= Yl U^,,^,{x)e'^'''+'^'^^^, 

(KoAd)eJ' 

where J is a finite subset o/Z x M and the coefficients f^reo,re£j ^^6 ^'^ H^(Kj^ ). 

The proof is given in section HI Observe that lA is periodic in ^o but almost periodic in (^0;Cd)- 
Proposition II .281 is applied to solve the equation 



£o(5eo, 5e,H' =[-{!- E) {L{d)Vl + D{0)Vp + (5.Z?(0)V°)V0)] 



e^(eo,?d) 



With this choice of Up we show in section [J] that the new approximate solution ^^(xj^o) in (|l-45p 
satisfies instead of (I1.43P the singular system 

(a) h,{U,) = O {e{K5 + C{5)e)) , 
(1.48) (6) il>{eUe)Ue - eG{x',eo) = 0{e^C{6)) on xa = 0, 

(c) Ue = Omt<0, 

where the errors in ()1.48p (a).(b) are measured in appropriate norms. Now one can apply our basic 
estimate (|2.55p for the linearized singular problem to conclude that the difference between exact 
and approximate solutions of the semilinear singular system (I1.18P satisfies for some constants C{5) 
and K: 

\Us{x,6o) -Ue{x,6o)\E= < K6 + C{6)e, for some s > — — . 

This estimate clearly implies the conclusion of Theorem 11.141 by choosing first 6 > small enough 
and then letting e tend to zero. 

1.6 Remarks on quasilinear problems 

In this article, we are able to rigorously justify a weakly nonlinear regime with amplification for 
semilinear hyperbolic initial boundary value problems. Our assumptions only deal with the princi- 
pal part of the operators, meaning that we only assume a weak stability property for the problem 
{L{d), B) obtained by linearizing at the origin and dropping the zero order term in the hyperbolic 
system. The weak stability is of WR type in the terminology of [BGR SZ02J . Despite the weak 
regime that we consider (O(e^) source term at the boundary and 0{e) solution), the leading profile 
equation displays some nonlinear features. We emphasize that the regime that we consider here is 
exactly one power of e weaker than the weakly nonlinear regime for the semilinear Cauchy prob- 
lem or for semilinear uniformly stable boundary value problems. As in [CGlOj . this power of e 
corresponds exacty to the loss of one derivative in the energy estimates. 

We believe that the techniques developed here can be extended to give a rigorous justification 
of weakly nonlinear geometric optics with amplification for quasilinear hyperbolic initial boundary 
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value problems of the form 

d 

dtv + Y,Bj{v)d,v + fo{v)=0, 
i=i 

(f,(v) = e^G (x', ^— ^"j on xa = 0, 
v = and G = in t < 0. 

The corresponding solution Vg; would be of amplitude O(e^). In particular the arguments used in 
Section [2] to obtain uniform estimates with a loss of one singular derivative for the singular initial 
boundary value problem might be extended to the corresponding singular quasilinear problem. 
There are however several new obstacles along the way, one of which is to extend the singular 
pseudodifferential calculus of [CGW12] in order to obtain a two-terms expansion of (1, 0) and (0, 1) 
compositions. The weaker scaling (e^ in place of e) should be sufficient to obtain the appropriate 
results. Let us observe that for O(e^) solutions, the principal part of the hyperbolic operator 
has coefficients that are uniformly bounded in W'^'^, which is precisely the regularity needed 
in [ Cou04t ICouOS ] to obtain a priori estimates and well-posedness. The leading profile equation 
obtained in this quasilinear framework is very similar to the one we have derived here, and we thus 
believe that a weak well-posedness result using Nash-Moser iteration should prove the existence 
of the leading proffie. For all the above reasons, we thus believe that the e^ source term on the 
boundary is the relevant "weakly nonlinear regime with amplification" in the quasilinear case, and 
we postpone the verification of the many technical details to a future work. 
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2 Exact oscillatory solutions on a fixed time interval 

2.1 The basic estimate for the Unearized singular system 

The goal of this section is to prove Proposition 12.21 below and its time- localized version, that is, 
Proposition 12. 9i These propositions provide the a priori estimates for the linearized singular system 
that form the basis for the Nash-Moser iteration of section 15.21 and the error analysis of section SI 
We begin by gathering some of the notation for spaces and norms that is needed below. 

Notations 2.1. Here we take s G N = {0, 1, 2, . . . }. 

(a) Letn :=mJ^\t\ Qt := Jln{-oo <t<T}, bn := M'^ x T\ b^r :=6fin{-oo <t<T}, 

and set ujt '■= ^+ H {— oo < t < T}. 

(h) Let H" = H^{bi^), the standard Soholev space with norm (y(a;', 0o))s- For 7 > 1 we set 
H^ := eT* H' and {V)s,^ := (e~^* V).,. 

(c) L^H" = L^(R+,H''{bn)) with norm \U{x,eo)\L2H- = |t/|o,s given by 

2 / \TT^' „ n \\2 



\U\o,s= \U{x',Xd,eo)\Hs^hn)dxd- 

The corresponding norm on Lp'Hi^ is denoted |T^|o,s,7- 

(d) CH^ = C{Rj^^H^{b^)) denotes the space of continuous bounded functions of Xd with values 
in H''{b9), with norm \U{x,9q)\ch- = \U\oo,s ■= sup^^yQ\U{.,Xd, ■)\H={bnT) ("^ote that CH^ C 
L'^H^). The corresponding norm on CH^ is denoted |T/|oo,s,7- 

(e) Let Mq := 3d + 5 and define C^^^^o •= (7(1^+, C*^°(6r2)) as the space of continuous bounded 
functions ofxd with values in C^^°{bQ), with norm \U{x, ^o)lco.Afo := \U\^cx,yyMQ,oo . Here L°°W'^^°'°° 
denotes the space L'^{M.+ ;W^''°'°°{bn))^ 

(f)The corresponding spaces on VLt «'"e denoted L^H^, L'^H^j., CH^, CH^j, and Cj^ ° with 
norms |t/|o,s,T> |f^|o,s,7,T; \U\oo,s,T, |f^|oo,s,7,T; o-nd \U\ o,Mo respectively. Onb^T w^ use the spaces 

H^ and H^j, with norms (C^)s,T CLnd {U)s,'^,t- 

(g) All constants appearing in the estimates below are independent of e, 7, and T unless such 
dependence is explicitly noted. 



The linearization of the singular problem (ll.lSp at U{x,9q) has the form 
(2.1) 



(a) ddUe + A{d,^ + ^]Ue+ V{eU) U^ = f{x, 9o) on ft , 



{b)B{eU)U,\,,=o = g{x',eo), 
(c) f/e = in t < 0, 

where the matrices B{eU), V{eU) are defined in (ll.ll|l^1 . Instead of (|2.ip . consider the equivalent 
problem satisfied by U"^ := e~'^^U: 

ddW + A ({dt + 7, 5.") + ^) il^ + V{eU) W = r(x, Oo) , 

f/'^ = in t < . 

^^The size of Mq is determined by the requirements of the singular calculus described in Appendix [XI 
^^Here and below we often suppress the subscript e on U . 
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Below we let A^) denote the singular Fourier multiplier (see ()A.2p ) associated to the symbol 



(2.3) A(X,7):= 7' + 



r + 



k(3 



e 



2' 



1/2 



X := e' + 



k/3 



£ 



The basic estimate for the linearized singular problem (j2.2p is given in the next Proposition. Observe 
that the estimate (12. 4p exhibits a loss of one "singular derivative" A^ . This is quite a high price to 
pay, which counts as a factor 1/e. In view of |CG10|. Theorem 4.1], there is strong evidence that 
the loss below is optimal. 

Proposition 2.2 (Main L^ linear estimate). We make the structural assumptions of Theorem \1.14\ 
and recall Mq = 3(i + 5. Fix K > Q and suppose \£ ddU\(jo,Mo-i + |C/|(^o,mo < K for e G (0, 1]. There 
exist positive constants £q{K) > 0, C{K) > and 7n (K) > 1 such that sufficiently smooth solutions 
U of the linearized singular problem (12. ip satisfur^l 

(2.4) iLT-^loo + ^^ < C{K) ^l^^/^IO'O + |V-'/''lo>o , {^D9"')o + i^x'g^h 



/or7>7o(K), 0<e<eo(i^)- 

r/ie same estimate holds if B{eU^ in ()2.ip is replaced by B{£U, eU) and V[eU) is replaced by 
V{£U,£U) as long as \£dd{U,U)\(jo,Mo-i + |f7,W|f;o,Mo < K for e G (0, 1]. 



Corollary 2.3 (Main Hl^^ linear estimate). Under the same assumptions as in Proposition {2^ 
smooth enough solutions U of the linearized singular problem (|2.ip satisfy: 

(2.5) lU^U, + m.,, + (^ < C(K) ( |A./--l.a + |V.-rl.a ^ (Wh4_(V^sZh\ 

\/7 V 7^ y/i y 

/or7>7o(K), 0<e<eo(^). 

Short guide to the proof. The proof of Proposition 12.21 is completed using the next three 
propositions, each of which has the same hypotheses as Proposition 12.21 In the first step of the 
proof of Proposition 12.21 we choose a partition of unity defined by frequency cutoffs XiiOj ^ = 
1, . . . , A^i + A^2, such that for z = 1, . . . , A''! the function Xi is supported near a point of the "bad" 
set T, while for i > Ni the function Xi is supported away from T. The estimates of Xi.D U^ for 
i > Ni are done in Proposition 12.81 For such indices, Kreiss symmetrizers in the singular calculus 
are used to estimate Xi,DU'^ without loss. 

To handle the piece Xi,DU^ ■, i < A^i> where the loss occurs, we also choose an auxiliary cutoff 
X^ in the extended calculus with the properties ()A.6p and use the decomposition 

(2.6) Xi,DU^ = Xh Xi,DU^ + (1 - Xh) XldU^ ■ 

The first term on the right in ()2.6p is estimated in the Proposition 12.41 while the second term in 
Proposition 12.61 Proposition 12.61 is the most difficult part in the analysis because in this frequency 
domain, some of the commutators are poorly controlled. 



^*Note that the norms \u\o.i and |Adu|o,o are not equivalent. 
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Proof of Proposition \2.2i . I) Partition of unity. The compactness of T (see Assumption 11.11]) 
and S allows us to choose a finite open covering of S, C = {Vi}i=i^,,,^Ni+N2 such that {Vi}i=i^,,,^Ni 
covers T and such that Ujy^^^^Vj is disjoint from a neighborhood of T. Since T C "H we can 
arrange so that for each i G {1, . . . , A'^i} there is a conjugator Qo,i(Cllj ^^'i diagonal matrix Bi^j(^) 
satisfying (jl.Sp in Vj. Moreover, we can choose a neighborhood O of (0, 0) E M^ such that for each 
i < Ni there are functions ai, Pi^i, and Pj^2 on O x Vj with the properties described in Assumption 
11.111 For these symbols, we shall use the substitution (t;i,i;2) — >■ sU{x,0o) to prescribe the space 
dependence. 

We let Xi(C)) ^ = I, ■ ■ ■ , Ni + N2 be a smooth partition of unity subordinate to C, and extend 
the Xi to all ^ as functions homogeneous of degree zero. We smoothly extend each Qo,j (as a matrix 
with bounded inverse) first to S, and then to all ^ as a function homogenous of degree zero. We 
take similar extensions in C of Pi^i, -Pi, 2) ^i,i and ai, but with homogeneity of degree 1 in the cases 
of Dij and cjj. As with Qo,j) the extensions of -Pj,i and Pi2 are taken to have bounded inverses rn 
Of course, for a given i < Ni the property (J1.12p is satisfied only for C/lCl ^ H- 

II) First estimate near the bad set. The first estimate deals with a piece of If^ that is 
microlocalized near the bad set T and in the frequency domain |(C',7)| ^ \kl3\/£. 

Proposition 2.4. Fix i such that 1 < i < Ni and let U^ := Xd Xi,DU"' where Xi ^-s defined in step 
I) above and x% is a cut-off in the extended calculus satisfying (|A.6p . Then we have the a priori 
estimate 

(2-7) \Ui lo,o + ^^ < C{K) \^ -2 + ^3/2 + — ^ + ^^j , 

/or 7 > jo{K). 

Proof of Proposition \2.4\ The loss of derivatives in the estimate prevents us from treating the zero 
order term P(eC/) Lf^ as a forcing term, as we would in a uniformly stable problem. Thus, we need 
to use an argument that simultaneously diagonalizes A and the lower order term T>{eU). 

We now set x'^Xi = X-: '^ '■= Xd^^ and look for an estimate v. We let A(X, 7) = —A{X,^) 
denote the singular symbol such that 

An=A({dt + j,d,.) + ^Y 

Dropping superscripts 7, we see from (12. 2|) that v satisfies 

ddv + Adv + V{eU) v = XDf + meU),XD]U = Xnf + r-i,DU , 
(z.oj 

B{eU) v\^^=o = Xd9 + [S{^U),XDp\xa=o = Xd9 + r-i,DU\xa=o- 

Here and below r_i,£) denotes a singular operator of order —1 (which can change from one occur- 
rence to the next) computed using the singular calculus. Similarly, ro,D will denote an operator of 
order 0. In spite of the loss of the factor Ad in the estimate (j2.4p . we are able to treat r-i^uU as 
a forcing term (see, for example, (12.17P below). A term like ro^^U/^ would be too large to absorb. 



Recall the notation C ~ {''' ~ ^Ti'?)- Sometimes we will also write (^ = (^',7) to match the notation of [CGW12] . 



^"Taking such extensions will reduce the number of cutoff functions we need later. 
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1. Simultaneous diagonalization. This diagonalization argument is similar to the one 
in |Cou04j . Let Qo{C) := Qo,i{C) and Bi(C) := ^i,i{C) be the matrices as in pTS]) such that 
(5o(C)''^(C)Qo'^(C) = lI5i(C) ill the conical extension of Vj. We define 

(2.9) w:=Qdv, 

where Q = Qo{X,j) + Q^i{eU,X,^). Here the matrix Q-i{eU,C,) is a symbol of order —1 defined 
for all C, but chosen so that on the conical extension of Vj, the matrix 

(2.10) Bo(ef/,C) := [Q-iQ^\lii] + QoV{eU)Q^^ 

is block diagonal, necessarily of order 0, with blocks of the same dimensions ni, . . . ,nj as those of 
Bi. Since the eigenvalues associated to the blocks of Bi are mutually distinct, a direct computation 
shows that Q-iQq , and thus Q-i, can be chosen so that the commutator cancels the off-diagonal 
blocks of QoT){sU)Qq . (The diagonal blocks of the commutator are all zero blocks and can 
therefore not cancel those of Qo V[eU) Qq .) Since QqK = I!)i Qq on Vj, (|2.10p implies the relation 

(2.11) QA + Qo^ = Bi<3+ [<9-l'^o^ID'l]<9o + <3o^ = BlQ + BoQo• 
Remark 2.5. [Entries of the symbol Q-i] Observe that the scalar entries of the matrix Q-i^d can 
be chosen to have the form 

(2.12) (Q_i,B)M = c(et/)a-i,D, 

where a_i(^) is of order —1 and independent of (x, 6), thus giving rise to a Fourier multiplier. 

Noting that x-dependence is absent in A and Qo and using the commutation property (|2.1ip . 
we have 

ddw = Qd ddv + {ddQ^i)DV = -Qd (A + 'D{eU))DV + Qd XdI + r^i,DU 
= -(Q A + QoV{eU))DV + Qd XdI + r-i,DU 
= -{BiQ + BoQo)DV + QDXDf + r-i,DU 
= -(Bi + Bo)dw + Qd XdI + r-i,DU . 

Here we have used the support property of x*^ to conclude (see Proposition lA.Qf a) and (b), as well 
as Remark lA.lOp 

(a) {ddQ-i)DV = r^iDU and 
(2.14) 

(6) (Bi Q^i)dv = Bi,B (Q-i)Df + r^i,DU . 

The fact that the statements in (|2.14p are no longer true when x'^ is replaced by (1 — x*^) i^i the 
definition of v is the main reason for the difference between Propositions 12.41 and 12.61 (the latter 
being much more technical). Here we are able to treat all commutators as forcing terms since they 
are order —1 operators. 

2. Outgoing modes. Recall that — Bi and — Bq are block diagonal: 
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so the system (|2.13p satisfied hy w = {wi, . . . ,wj) can be written as a collection of J decoupled 
transport equations 

(2.15) ddWj = {iujj)DWj + Cj^DWj + ro^of + r^i^oU 

with Ke{i(jjj) < —07 for 1 < j < J', and Re{iujj) > 07 for J' + 1 < j < J (c > denotes a 
constant). 

Following the strategy of |Cou04] , we now give two preliminary estimates of the outgoing modes 
Wj, j > J' + 1. Taking the real part of the L^(il) inner product of (j2.15p with A'jjWj, we obtain 

^^~2 ~ ^ ^^ ^^^ i^^j)D'Wj,ADWj)L2(^Q) + Re {Ad Cj^dWj,AdWj)l2^q) 

+ Re (Ad Qd Xofi ^dWj)l2(^^) + Re (Ad r^i^^U, AdWj)l2(^^^ . 
Since Re (iujj) > 07, we get after absorbing some terms on the left 

|2 , /A _„.. /n\\2 ^ ^ f\\ _ r\2 , irV|2 



(2.16) 7 \^DWj\lo + i^DWjmo < - [\^Df\io + |f/|0,0 
Here, for example, we have used Young's inequality and estimated 

(2.17) |Re(ADr_i,D^,ADu;j)i2(o)| < —\U\lo + 6-f\ADWj\lo , 

with (5 > small enough so as to absorb the term lADit'jIoo from right to left. Similarly, taking 
the real part of the L^ inner product of (|2.15p with Wj on [x^, 00) x b^ instead of 0, we obtain for 
all Xd > 0: 

(2.18) 7 kjlo.o + {wA^d))l < ^ (l/lo.o + ^ |f>lo,o) • 

3. Incoming modes. Estimating Wj for j < J' in a similar way, but now using Re {itOj) < —07 
and pairing the corresponding transport equation with Wj, we obtain 

l\w,\lo<C{w,ml + ^ {\f\lo + \r-i,Dif\lo)<C{w,ml + ^ ^\f\l^ + l^\U\l^^ . 

Taking the real part of the L^ inner product of ()2.15p with Wj on [0, Xd] x bQ rather than on all fi, 
we also obtain 

(2.19) ^|^^.|2^^ + (^.(3;,))2<c(u;,(0))2 + ^ (^|/|2^^ + J_|{7|2^^^ . 

4. Boundary estimate. We observe that v can be expressed in terms of w as 

(2.20) v = {QQ^)DW + r.i^DU- 

Recalling the boundary condition in ()2.8p and using the decomposition Qq {() = [Qm(C) Qout{0]i 
we let accordingly w = {win-,Wout) and rewrite the boundary condition in (j2.8p as 

(2.21) B{£U) Qin,DWin\xi=Q = -B{eU) Qout,DWout\xa=0 + Xd9 + r-i^DU\xa=0 ■ 
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By prT2|) we have on H 

B{eU) Q,n = P,-' {Pi B{eU) Q,„ P2) P^' = Pf ^ (^'' ^^ + '""^ ^ P^', 
so using the rules of singular calculus, we get 

Ai:)B{eU)Qin,DWin\xa=0 = (.Pi^)d i ^ ^ jj {P2^)DWin\xa=0 + ro,DWin\xa=0- 

With dOH), this implies 
(2.22) 

{iPl')D (^ ^ ''''' ^^ j] iP2')DW^n\x,=0)0 < C [{ADWout\x,=o)o + (Ad5)o + (f>U,=o)o) • 

We have Pi^d{Pi )d = I + t-i,d so up to choosing 7 large (and absorbing the r_i_£) term), the 
estimate (j2.22p implies 

(2.23) ( (^ + '""^ ^^ \ {P2^)DW^n\x,=0)Q < C ((AD«^oniU,=0)0 + (Ad5>0 + (^U,=o)o) • 

Letting ( / ) := {P2 )DWin\xi=o we find, using the fact that a is real and choosing again 7 large 



enough. 



^ ( ^ ^ "" A. /) b) >o ^ h (^' ^-^^0 + ^^--'^o) > ^ (-1, -')g. 



Thus, from (j2.23p we may conclude 

(2.24) 7 (t(;jnUd=o)o < C [{ADWout\xi=Q)o + (Ad5')o + {U\xi=o)oj ■ 

5. Conclusion. Combining the estimates (f2l^ . ([2J8]) . (fZT9]) . and (jTM]) we obtain 

7^ \Win\lfi + 7^ kinlL.O + 7 |AD?yo«t|o,0 + (ADU^OMtU'd=o)o + 7^ ko«t|L,0 

< C7 (^^^ + (A,,)g + Mm + (^|,^^„)g^ . 

Recalling the relation ()2.20p between v and w = Qdv, we obtain the estimate (12. 7p for 7 large 
enough. Though we have not emphasized the regularity of the symbols, we claim that the regularity 
assumption on U is sufficient to apply the rules of singular symbolic calculus at each step of the 
above calculations, see Appendix |Al We feel free to skip some of the details. D 

III) Second estimate near the bad set. Next we estimate the second term on the right in 
(12. 6j) . In the statement below Qd = Qo,d + Q-i,D is the same operator as that constructed above 
in the proof of Proposition 12. 4i Since {Qo,d)~^Q-i.d has norm less than one as an operator on L^ 
for 7 large, we can define {Qd)~^ as an operator on L^ using a Neumann series. 
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Proposition 2.6. Fix i such that 1 < i < Ni, let C/g •= (1 ~ X%)Xi,DU'^ (md write 
(2-25) U^, = Ul^ + Ul^^,, 

(2.26) f72% := (Qi^)-^K„,0) and Ul^^, := iQDr\0,Wout). 

Then we have 

(2-27) \U^^Jo,o + ^ + + ^575 

^1^ .-i,V7 I , (V:.'A^^Lr,^^,,U,=o)o ^ |(v,,,^)/7|o^o ((V.^7)g^)o 



D ^2,mlU,U -r ^ - w ^ ^2 ^ ^3/2 

72 + ^3/2 ) • 

Proof of Proposition [KR 1. Simultaneous diagonalization ^vith a new remainder. We now 

set (1 — x'^)Xi = X) ^ := Xd^^ — ^2^ ^ •= Qdv = {win,Wout) and repeat the proof of Proposition 
12.41 down to line (|2.13|) with only one change. Suppressing superscripts 7, we now have instead of 
(I2J3I1 

(2.28) ddw = -(Di + Oo)dw + ro^of + t-i^dU + R'hv, 

where i?^ = i?|) — i?^ with operators R^^^R'fj defined in (j2.30p below. As can be checked by 
looking at each line in ()2.13p . we also observe that the remainder term r^i^£,U reads 

(2.29) Qd [V{eU),XD\U- iQo,D Bo,d - (Qo Oo)d) XdU - Q-i,d V{eU) XdU + ^o,d Q-i,d XdU ■ 
The equalities of ()2.14p should be replaced by 

^2 3Q^ («) {ddQ-i)DV = RdV , and 

ib)ilD)iQ^i)DV = Bi^D{Q-i)DV + Rhv, 

where in view of Remark 12.51 the scalar entries of i?^ and i?^ have the forms 

(2.31) {ddc{eU)) a-i^D and [qi,d, c(e[/)] a_i,z) , 

respectively. In (I2.3ip . ai(C) denotes one of the diagonal entries of Bi(C). Here and below O-i^d 
denotes a singular operator of order —1 associated to a symbol a_i(C) which may change from term 
to term. We claim 

C C 

(2.32) \Rdv\o,o < - blo,o and \V ^> Rf)v\o,o < C |V^/A^^v|o,o + - |v|o,o- 

7 7 

Indeed, the first estimate is clear from the definition ()2.3ip since ddc{eU) is a bounded function 
and [ai_/),c(e[/)] is bounded on L? by the rules of singular calculus. The L^ estimate of Vx'R\)V 



'Here Win £ C and Wout £ C^ ^ are defined by the same diagonalization procedure as in the previous proof. 
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is readily obtained by differentiating the product {ddc{eU)) a^i^^v with respect to x' so we now 
focus on the estimate of V^'R^v. Setting vq^d := [oi^d, c{eU)], which is a bounded operator on L^ 
by the rules of singular calculus, we have 



l^^x' Rhvlofi < \ro,D^x' a-i,Dv\o,o + \[^x',ro,D]a-i,Dv\o,o < C jV^^-'A^, u|o,o + 



7 



where we have used Proposition IA.9l fc) to estimate the double commutator term [V^^/jTo^d]- 

The estimates (j2.32p will be used below to handle the source term Rj^v appearing on the right 

of (|2.28|) . Unlike what happened in the proof of Proposition 12. 4^ (|2.32p seems to be the best we 

can hope for the bad (1, —1) product. 

2. Incoming modes I. In place of (I2.15P we now have 



(2.34) 7^ \win\lo + 7' km|L,0 < 7' {mn\x,=0)l + - (|7 /lo,0 + l^lo,o) 



(2.33) ddWj = {iujj)DWj + Cj^dWj + ro^of + t-i^dU + R'h'^, j = 1, • • • , J 

and for the incoming part Win we obtain as in ()2.19p : 

C 

7 

This L^ estimate does not cause any problem because we have the good L^ control (|2.32p of the 
additional term R'^v appearing on the right of (|2.33p . 

We can also repeat the same argument as in step 4 of the previous proof for the boundary terms 
because these estimates did not rely on the support properties of x'^ (the support properties of x^ 
were only used to control some error terms in the interior equation). We thus derive again 

7 (w^mUd=o)o < C {^{ADWout\xd=o)o + {^D Xd9)o + {U\xd=o)o) ■ 

Using the fact that on the support of 1 — x^(C') ^ , 7), we have 

i^,7i<cie',7i, 

we obtain the estimate 

{ADWout\xa=o)o + {AdXd9)o < C {{Vx' ,^) Wout\xa=o)o + C {(V^' ,j) g)o ■ 
At this stage, we have thus derived the bound 

(2.35) 7' km|g,o + 7' Iwinllfl < ^ (l7 /lo,o + |f>lo,o) + C (((V,,,7) ff>o + {U\x,=o)l) 

+ C ((Va;', 7) Wout\xa=o)o ■ 

3. Outgoing modes. Let j £ {J' + 1, . . . , J}. Taking the real part of the L^(Q) inner product 
of (|2.33p with —d^^Wj for /c = 0, . . . , d — 1, and with 'y'^Wj, we obtain in place of (|2.18p 



C 
7 
Here in place of (|2.17p we have used the estimate 



(2.36) i\(yx',l)wout\lo + \(yx',l)wout\'L,o<- (\i'^x',l)f\lo + \U\lo + \dx,A^^U\lo) ■ 



\{d:rk r-i,DU, d^^Wj)L'2^n)\ < -^ \dcc^ A^^f7|g_o + — l^lo,o + ^lldx^Wjl^Q, 



7 7 
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along with a similar estimate for \{dx^R'£)V,dxf,Wj)£^2rQ\\. Such estimates follow either from the 
precise expression of the remainder dx,^ r^i^uU (differentiating each of the four terms in (j2.29p with 
respect to Xk) or from ()2.32p . As usual, we then choose S small enough so as to absorb terms from 
right to left. 

Adding (J2.35P for the incoming modes and (|2.36p for the outgoing modes, we have already 
shown how to control the first line in (12.270 by the quantity on the right of the inequality ()2.27p . 
At this stage, it thus only remains to control the term V^/ A^ ^2in- 

Remark 2.7. At this point we can see the need to estimate the remaining terms on the left in the 
estimate (j2.27p as well as the similar terms on the left in the Kreiss estimate (j2.43p below. We must 
estimate those terms in order to be able to absorb the terms involving V^;' A^ U^ on the right side 
of (j2.27p . We recall that these terms come from the "bad" composition term Bi^/j (5-i,d and from 
ddQ-i,D- Observe that there is no need for a separate estimate of Vx' A^ U^ since, for example, 

|V,.A^iL^7|o,o<C|f77|o,o 

because \X,^\ > C \C',j\ on the support of x^- 

4. Incoming modes II. Here we begin to estimate the terms in the second line of (|2.27p . 
The problem satisfied hj v = U2 has the form (]2.8p with the new cut-off function X = (1 — x*^) Xi- 
We now introduce the functions v := A^ v and v' := dx'V, where dx' denotes any of the tangential 
derivatives dxQ, ■ ■ ■ ,dx^_^, and see that the function v' satisfies 

ddv' + Adv' + Vi£U)v' = dx> A-1 xnf + dx' A^' [V{eU),XD]U 
(2.37) + dx'{{V{eU) - A'j^)'^ V{eU) Ad)^) - dV{eU) ■ edx'Uv , 

B{eU) v'\x,=^ = dx' A^i Xd9 + dx' A^' [B{eU), Xd\UU,=g + dx' [S(eC/), A^>U,=o. 

We can thus diagonalize the problem for v' with the same operator Q£, as before. Introducing the 
function w' := Qdv' , we find that w' satisfies 

(2.38) 

(a) ddw' = -(Pi + Bo)dw + Qd dx' A'^ Xnf + - ro,DV + - r^^oU + Qd dx' A^^ [V{eU),XD\U , 

-^ 'y 'y 

{h) B{eU) Qin,Dw[^ = -B{eU) Qout,Dw'^t + dx' A^' Xd9 + dx' A^' [B{eU),XD]U 

+ dx'A-^^[AD,B{eU)]v + r^i,DV. 

where we have collected many terms (for instance the i?^ operator as in (I2.28P ) into remainders of 
the form 7"^ rQ^D. For instance, we have used 

Qd dx' {{V{eU) - A^^ V{£U) Ad)v) 

= Qd {V{eU) - A^^ V{eU) Ad)v + Qd {dx'V{eU) - A^' dx'V{eU) Ad)v 

= - ro,DV + - roDV = - tqdv + - tq^dU . 

7 7 7 1 ' 

Next we fix an index j G {!,... ,</'}• Taking the real part of the L^(il) inner product of 
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(j2.38p (a) with w',, we obtain the standard L^ estimate for incoming modes: 

+ ^\d^,A-^'[VieU),XD]U\lo- 

The last term on the right of the inequahty is estimated by similar techniques as above, namely we 
write 

d,. A^i [VieU),XD] = meU),XD] d,' A^' + [d.' ^~D\V{eU)] Xd " Xd [d.' ^n'M^U)] , 

and then decompose the commutator [d^' A^ ,T){eU)] as follows: 

[a,. A^\V{eU)] = d^> (A^i V{eU) Ad - V{eU)) A^' + {d^,>V{eU)) A^^ . 

In the end we obtain the estimate 

(2.39) ^\u,'Jl^<C{w',^\,,=o)l + ^ {\d.'f\lo + \d.'AD'U\lo + \U\lo) > 



and we thus wish to control the trace of w'^^. 

5. Control of the trace of w'^^. Using (|2.38p (6). and arguing as in step 4 of the proof of 
Proposition 12.41 we obtain the boundary estimate 



(2.40) 7 {w',J,,=o)o < C [{Adw',J^,=o)o + {d,,g)o + {d,, [B{eU), Xd]U\.,=o)o^ 

+ C {{d^'[An,B{eU)]v\,,=o)o + {v'U,=o)o) . 

The two commutators appearing on the right of (j2.40p are dealt with as in the previous step and 
we get 

7 {wL\xa=o)o < C ((ADW)ontUrf=0>0 + {dx'9)o + {U\xa=o)o + {dx' ^D^U\xa=o)o) ■ 
Combining with (|2.39p . we have derived 

(2.41) 7 \w'Jl, + {w'J,,=o)l < ^ {AdKJx,=o)1 + ^ {\dx'f\lo + 1^:.' ^~dU\Io + l^lo,o 

C 



+ - {{dx'gfo + {U\x,=o)l + {dx' ^d'U\x,=o)1 



r 

We expect Ajjw'^^^ to be comparable to d^'Wout and thus use (j2.36p : this is checked and made 
precise in the next step. 

6. Relation between Aoui' and d^'W, and conclusion. Using the definitions 

w' = Qd v' = Qd dx' A^ V and w = Qdv , 

and the fact that A/j commutes with Qo,D: we compute 

Adw' = Qd dx'V + ro,Dv' = d^'W + tq^dv + vq^dv' ■ 
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We have thus derived the bound from above 
1 C / 

which we combine with ()2.4ip and (j2.36p to obtain 

(2.42) 7 \w'Jlfl + (*LU.=o)g < ^ (|(9x', 7) /lo,o + l^.' A^'f^lo.o + |f>lo,o 






+ -2 {i^x'gfo + {U\.,=o)l + id.' AD'f>U.=o)o) 



It only remains to derive a bound from below to go from w^,^ to dx'A^ ^2 in- ^^ ^'^s.i observe 
that estimating dx'A.]^ ^2^in ^^ claimed in (j2.27p amounts to estimating Qndx'^J) ^2in- ^^ ^^^ 
the relation 

and we use the special "decoupled" form of the coefficients of Q-i to show that the commutator 
[5i.'A^\(5_i,£)] reads 

[dx'A]j^,Q-i^D] = — ro,D + — ro,D {dx' A^'^) . 



Similary, we can write 



so we obtain 



w[n = dx'Aj^Win + — tq^dv + — ro,D {dx' Aj^) 



Qd dx'Aj]-Ul^^ = w[^ + — ro^^il + — ro,D {dx' A^^)U . 



We have therefore proved that (j2.42p implies that the second line in (|2.27p is controlled by the 
terms on the right of (j2.27p . This completes the proof of Proposition 12.61 D 

IV) Estimate a'way from the bad set. The next proposition provides a Kreiss-type estimate 

for the terms Xi,DU''', where i > Ni. 

Proposition 2.8. Fix i such that Ni + I < i < N2 and let (l^ := Xi,DU'^ ■ We have 

(')A'i\ ITT^I I (^3Ud=o)o , IV7 A-lry^l I i^x'Ajj U^\x^^q)q 
{2A6) (^3 0,0 H -;= h Vi:/ Arj L/3 0,0 H ^ 

^/7 Vl 

^ ^ f \r\o,o + \^x'A^'P\o,o ^ {g^)o + (V,.. A^'g^)o 
~ \ 7 \/7 

, \U^\o,o + |V,. A^'iPjofl ^ {U^\x,=o)o + jVx' A^'U^\x,=o)o \ 



7^ 
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Proof. 1. L^ estimate. The first step is to prove the Kreiss-type estimate 



\/7 ~ V 7 \/7 7^ 



(2.44) |L^3^|o,o + ' ' 'T < C 1 ^4^ + ^ + ^4^ + 



For this we define the good set G C S to be a neighborhood of the closure of U^J^ +i^i such that 
G is disjoint from T; here the uniform Lopatinskii condition is satisfied. The classical construction 
of Kreiss symmetrizers [KreTOt rCP82| provides us with an A^ x A^ symbol R{C), homogeneous of 
degree 0, such that for some positive constants C, c and C/lCl G G we have 

(a) R{C) = RiCr 

(2.45) (6) -Re(i?(C)A(C))>c7/iv 

(c) R{O + C13{0)*B{0)>cIn. 

We take a smooth extension of R to all C as a symbol of order such that p.45p (a) holds. Observe 
that by continuity (j2.45p (c) implies 

(2.46) R{C) + CBieUy l3{eU) >cIn for e small enough. 

As observed in [Wil02j we may now use Rd, the singular Fourier multiplier associated to the 
symbol R{X,^) as a Kreiss symmetrizer for the singular problem. Let Xi = X^ v •= Xd^^'i ^^'^ 
denote by (•,•) the L^ inner product on b^. Using the equation ()2.8p to expand dd{v,RDv) and 
integrating in x^ over [0, oo) we obtain 

- {v\xi=o,RDv\xi=Q) = -2Re {RdAdv,v) 

- 2 Re {Rd V{eU)v, ^) + 2 Re [Rd XdP, v) + 0{\U"'\loh'')- 
From (I2.45p (b). (I2.46P and the localized Carding inequality (Proposition lA. iTl) . 



(2.47) Re{{R + C B{eUy B{eU))DvU,=o, v\x,=o) > c {vU,=o)t - C ■ 



7 



we easily derive the estimate ()2.44p . 

2. Estimate of V^-zA^ IJ^. Set v := A^ v and for now let v' denote one of the tangential 
derivatives djV, j = 0, . . . ,d— 1. Then v' satisfies the system (j2.37p . where the truncation function 
X has changed but the forcing terms have exactly the same expression. An argument just like the 
one that gave the estimate ()2.44p yields 

I V^/ A^ U^ |o,o H ;= < G I ■ h 



\/7 \ 'y V^ 



l^lo.o + |Vx'A^^;7^|o,o (^^U,=o)o + {Vx' K-^'W\x,=o)o 

T 7 

Here instead of (|2.47p we have used 

Re m . CBWr BieUM.^^, .',,.0) , c (.'U,.„,§ - C ^^^.Ip,^ , 
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to recover the estimate of the trace of v' . The I? estimates of the forcing terms in the interior 
and on the boundary are exactly the same as in steps 4 and 5 of the previous proof. Since these 
estimates are actually simpler or similar to those we dealt with in Proposition 12.61 we feel free to 
skip the details that should become more or less familiar to the reader at this stage. D 



V) Conclusion. We use the previous propositions to complete the proof of Proposition 12.21 
It remains to estimate |f/'^|o,i and |f7'''|oo,o- Summing the estimates (|2.7p . (|2.27p . and (j2.43p 
over i G {1, . . . , A''2} and absorbing error terms from the right by taking 7 large, we derive 

(2.48) 



V7 Vt V 7^ 7^/^ 

where we have "forgotten" on the left of the inequality the additional control of V^./ A^ f/'^ (this 
term has played its role, meaning that it was used to absorb some bad terms appearing on the 
right). This gives exactly (j2.4p with the additional control of f/^ + IJl in U^(l?\ This additional 
property will be used in the proof of Corollary 12.31 below. D 

Proof of Corollary \2.kA We first estimate the first order tangential derivatives. We can apply the 
a priori estimate (12. 4p to the problem satisfied by dt^', 60)^^1 which is obtained by differentiating 
(12^ . This yields 

(2.49) \U^\o,i + — < C{K) ^ -2 + -^^ 

which is the same as (|2.5p . except for the absence of |f7'''|oo,o on the left. Here we were able to treat 
commutators as forcing terms because, for example, 

[V{eU), 9(,,,eo)] U^ = -(dP(eC/) • e9(,,,,,)[/) U^ , 

and the factor of e coming out from the commutation allows us to estimate 

\AD[V{eU),d^,,^eo)]U^\o,o < C\U^\o,i. 

It thus only remains to estimate the norm |{/'''|oo,o- For (52 > to be chosen, we take < 5i < 62 
and consider a symbol of order zero in the extended calculus, x'^i^'i )7)) such that 

< X' < 1 , 

X U '— '7 1 = 1 on <^ 1^ ,7| < (5i ^— 



suppx C <^ |^,7| < 02^— } 



We then write V^ = x\)U^ + (1 ~ x\))U^ and begin by estimating |(1 — x%)U^\Q,oa by using the 
Sobolev-type estimate 

(2.50) 
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Using the equation (|2.2p and the fact that 



we obtain 



1(1 - xh) ddii^Wo < |Ad (1 - xhWko + 1(1 - xh) ^u^ko + 1(1 - xh)P 



10,0 



<C |f7^|o,i + |r|o,o <C |f7^|o,i + 



7^ 



where the last inequality follows from |/'*'|o,o < C'|/'^|o,i/7. With (I2.50J) this gives 

|AD(r)io,i' 



(2.51) 



\il-xh)U^\oo,o<C{\U^\o,i + 



r 



To estimate \x%U^\oo,o we observe that since /3 G T, we have for ^2 > chosen small enough 



x^(e',^,7)=x^(e',^,7)E^«(^'^) 



for the Xi chosen in step I) of the proof of Proposition I2.2i Thus, 

W'hU^oofi < Ixh {U2 + ^2^)|oo,0 <m + f/2^|oo,0 , 

with U^ and U2 defined in Propositions 12.41 and 12. 6i We can then apply the a priori estimate ()2.48p 
and obtain 

I e,V7| ^^^ lAprlo.O + IVx'rlo.O , (ADg^)0 + (V.'g^)0 ^ 
\XDU'\oofi < C (^ -^^ + J . 

With (f23T]) and K^ . this completes the proof of Corollary EH D 

Let us quickly observe that the genuine Garding's inequality was used only once, in the proof of 
Proposition [2]2l namely in ()2.47p . In all other cases, we only used Plancherel's Theorem for Fourier 
multipliers. This explains the slight difference between ()2.43p and (12. 7p . ()2.27p for the powers of 7. 



Next we "localize the estimate" to Q,t- Sino 



IAd^ 



10,1 



i,d.' + ^^]r 



0,1 



7,a., + ^)/ 



0,1,7 



we can rewrite the a priori estimate (|2.5p for solutions to the linearized system (j2.ip as 
(2.52) |f7|^^o^^ + |[/|o^i_^ + 



{U\xa=o)l,-y 



^ 



m, 



< ^^^^ / 1(7,9.' + ^)/|o.i,7 + |Vx'/lo.i,7 ^ {il,d.' + ^)5)i,7 + (V.'5)i,7 



7' 



.3/2 



Here "~" denotes equivalence of norms with constants independent of e and 7. 
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Suppose now that the singular problem (j2.ip is posed on Qt instead of il. Given / G L^H^, 
one can define a Seeley extension / G Lp'H^ such that 



7> ^x' + 



Pde 



f 



+ |V,,/|o,i < C 



0,1 






+ |V,//|o,l,T , 



o,i,r 



where C is independent of 7, e, and T. It is readily checked that the same extension satisfies 



(2.53) 



7, dx' + 



/3< 



/ 



+ \^x'f\o,l,-y < C 



0,1,7 



7, (9^./ + / 



+ |Vi.//|o,l,7,T 



0,1,7,^ 



where again C is independent of 7, e, and T. We claim that changing /, g, and U in {t > T} 
does not affect the solution of (I2.ip in {t < T}. (This causality principle is discussed further below 
together with the existence of solutions to the linearized system (j2.ip .) Hence the estimates ()2.52p 
and (I2.53P imply the following estimate for the singular problem on il.T- 



\U\ 



oo,0,7,T 



+ |f^|0,l,7,T + 



(^Ud=o)i,7,r 



V7 



f^a, 



<C{K) 



1(7,5.' + ^)/|0,l,7,T + |V,,/|o,l,7,T , ((7,5x' + ^)5>1,7,T + (V,.ff)i,7,r 



r 



7 



3/2 



Let us now consider the linearized singular problem (j2.ip on Qt with data of the form ef, eg 
instead of / and g. We note that 



7, d^' + '- ef 



< C'|/|o,2,7,T and ( 7,5./ + 



0,1,7,T 



(3de, 



£g)l,-y,T < C{g)2,^,T- 



Let us write the linearized operators on the left sides of (12.1j) (a) and (b) as L,'{eU)U and M'{eU)U 
respectively, and define 

4([/){7 := - V{eU)U , B',{U)U := - M'{eU)U. 

We have proved: 

Proposition 2.9. Fix K > and suppose \eddU\ o,mo-i + \U\ o,mo < K for e G (0, 1]. There exist 
positive constants eo{K), 7o(-f^) such that solutions of the singular problem 

C'^{U)U = f onQr, 
(2.54) B'^{U)U = g on b^T , 

il = in t < 0, 
satisfy 



(2.55) 



|t^|oo,0,7,r + |^|o,l,7,T + 



(f^U'd=o)l,7,T . ^/j^A /^l/lo,2,7,T , (5)2,7,T 



<C{K) 



+ 



v3/2 



\/7 y ^^ ^o 

for < e < eQ{K), 7 > jq{K), and the constant C{K) only depends on K . 

The same estimate holds if B{eU) in (|2.ip is replaced by B{eU,eU) given in (|1.9p . and T>{eU) 
is replaced byV{eU,eU) given in (jl.lOp . as long as there holds \edd{U,hl)\ q,mq-\ + \U,hl\ o,mq < K 

/or EG (0,1]. 
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2.2 Well-posedness of the linearized singular equations 

In this short section, we explain why the analysis in |Cou05j gives existence and uniqueness of a 
solution to the linearized singular problem (|2.54p for which the estimate (j2.55p holds. First of all, 
we can define a dual problem for (j2.ip that reads 



, , dJj-rk*{d^.^^^^ij + V(eU)\J = f{x,eQ) onO, 

(2.56) V e / 

A^(ef/)f7U,=o = 5(a;',^o), 

where A* is obtained from A by first multiplying the system by the constant matrix 5^;, then by 
integrating by parts on il and eventually by multiplying by (-BJ)^^. The zero order term is also 
changed accordingly. Following the standard procedure described for instance in |BGS071 Chapter 
4.4], the matrix M giving the adjoint boundary conditions is chosen such that for all v sufficiently 
close to the origin, there holds 

where B\{v) and M.\{v) are additional matrices depending smoothly on v. 

The expression of A* shows that this singular operator coincides with the operator obtained by 
applying the substitution d^' — >• d^' + fiddf^ /e to the dual operator 

d 

i=i 

It is known from the analysis in |BGS07l Chapter 8.3] that the latter constant multiplicity hyper- 
bolic operator with boundary conditions given by M{v) gives rise to a boundary value problem 
in the "backward" WR class (one just has to replace 7 by —7 for this dual problem). When we 
apply the singular transformation d^' — >• d^' + (ido^/e to the boundary value problem defined by 
(-Lo(9)*, A^(eC/)), we can reproduce the analysis of the former section and show that the same type 
of a priori estimate as in Proposition 12.21 holds for (|2.56p . 

For all fixed e > small enough, we have thus proved that both the forward problem (j2.ip and 
its dual problem ()2.56p satisfy an a priori estimate with a loss of one tangential derivative. The 
estimates depend very badly on e because singular derivative d^' +(3dgg/e is estimated by 1/e times 
the tangential H^ norm with respect to {x',6q). Nevertheless, we can at this stage reproduce the 
arguments of |Cou 05] to show the existence and uniqueness of L^ solutions to (j2.ip when the source 
terms / and g satisfy f,dgfjf,dx'f G L^(Or), g S H^{b^T)- The analysis is actually much simpler 
than in |Cou05j because most of the technical difficulties in |Cou05] arise from commutations with 
the hyperbolic operator. Here the hyperbolic operator has constant coefficients so commutation 
with any scalar Fourier multiplier is exact. The analysis in |Cou05j also shows that weak solutions 
are limit of strong solutions when the hyperbolic operator has constant coefficientcfl so we can 
show that weak solutions satisfy the energy estimate ()2.4p with constants that are uniform with 
respect to the small parameter e. Such global in time estimates imply the causality principle that 
"future does not affect the past" and can be localized to i^T by the extension procedure previously 
described. 



^Weak solutions are only "semi-strong" solutions when the hyperbolic operator has variable coefficients. 
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2.3 Tame estimates 

In this section we prove higher derivative estimates for the hnearized singular problem (j2.1|) . 
first in the "pre-tame" form of Proposition 12.111 and then in the final "tame" form of Proposition 
12.161 which is suitable for Nash-Moser iteration. Propositions 12. 12] and [2rT5] give pre-tame and tame 
estimates for second derivatives. 

Notations 2.10. (a) Let L'^W^''^ = L°°(E.+ ,W'^'°°{bn)) with norm \U\loowi,oo := \U\* . We also 
write \U\Loo(^n) = \U\^, {V)L^(^bn) = {V)^, {V)wi,^^hn) = i^)* , \U\L^(nj,) = \U\^,t, etc.. 

(b) For /c G N, let d denote the collection of tangential operators d?, g -, with \a\ = k (a is a 
multi-index). Sometimes d^ is used to denote a particular member of this collection. Set d^<j) = (f). 

(c) For k £ {1, 2, 3, ... }, denote by d^ ' (f) the set of products of the form (9"^0jj) • • • (c?"''0i,,) 
where 1 < r < k, ai + ■ ■ ■ ttr = k, Qj > 1. Set d^^'cf) = 1. 

(d) For r > 0, let [r] denote the smallest integer greater than r. 

Our first goal is to prove the following "pre-tame" estimate for solutions to (|2.54p . 

Proposition 2.11. Fix K > and suppose \eddU\QO,Mo-i + \U\qo,mo < K for e G (0, 1]. For s > 
in any fixed finite interval there exist positive constants eQ{K), 70 (-f^) such that the solution to the 
linearized singular problem (I2.54P satisfies 

(O r^7\ IrVl I IrVl , {U\xa=o)s+l,-y,T 
[Z.bl) |t;|oo,s,7,T + |<^|0,s+l,7,T H 1= 

^ ^(^) (^ -2 + ^3/2 + -2 + ^57^ j ' 

/or < e < eQ{K) and 7 > 7o(i^)- 
Proof. The problem satisfied by 9*C/ is 

C'^{U)d'U = d'f + - [V{eU),d']U, 

B',{U)d'U = d'g + - [B{eU),d']U. 

In applying the estimate ()2.57p to this problem we must, for example, compute d'^{\p{eU),d^]U), 
which is a sum of terms of the formal 

V{eU) d^^^ {eU) d^U , where j + A: = s + 2, j > 1, 

and T) is some smooth function of its argument. Since j > 1, we can rewrite this as 

V{eU)d'^^~^\eU)d{eU)d''U. 

Using Moser estimates we obtain 



- V{eU) d^^-^^ {eU) d{eU) d^U 



< C{K) |t/|.,T \UWs+2n,T + C{K) |t/|o,.+l,7,T. 
0,7,T 



The contribution from the final term on the right can be absorbed by taking 7 large enough; thus 
this explains the third term on the right in (j2.57p . The final term on the right in ()2.57p arises by 
the same argument applied to the boundary commutator. D 



^''More precisely, each component is a sum of such terms. 
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Next we prove estimates for the second derivatives 

4'(C/)({7», tl") = d,V{eU) ({7^ il'') , 
B'^{U){U'', if'') = dyl3{eU) ({/", if'') , 

where we use the notation 

N 

d,V{eU) {U\U') := ^ {d,p{eU) U^) U\ 

i=l 

Proposition 2.12. We have 

{a)\C:{U){U\U'')\^^sa,T< 

C{\U\^,^T)\\U"'\oo,s,-y,T\U \*^T + \U |oo,s,7,T|f^'^|*,T + £|^|oo,s,7,T|f^'^|*,T|f^ |*,TJ , 

(6)|4'(^7)({7^f7^)|o,s+l,7,T< 

C(|C^|*,t) ( |f^°|o,s+l,7,T|f^ \*,T + \U |o,s+l,7,T|f^"|*,T + £|t^|o,s+l,7,T|f^'^|*,T|C^ \*,T 

{c){B'l{U){U\U'))s,,,T< 

Proof. For t < s one computes d^{C"{U){U"', f/^)), which is a sum of terms of the form 

V{eU) a<''> (eU) d^if" a"f7^ where k + l + m = t. 

Thus, the first estimate follows directly from Moser estimates. The remaining estimates are proved 
the same way. D 

In the iteration scheme of section 15.21 we will use H^ spaces on the boundary, while in the 
interior we use the following spaces. 

Definition 2.13. For s G {0, 1, 2, . . . } let 

E^ = CH^nL'^H^~^\ with the norm |C/(x,^o)b|, := \U\oo,s,T + \U\o,s+i,T 

Ef^^T = CH-;^T n L'^H^'jr\ with the norm |C/(x,eo)b^,^ := |C^|oo,s,7,T + \U\q,s+i,^,t- 

Remark 2.14. By Sobolev embedding we have 

s > [{d+l)/2] ^ E^C CH't C L°°{nT) 

s>[{d+ l)/2] + I^E'tC CH't C L°°(1+, P7^'°°(61]r)) 

s>[{d+ l)/2] + Afo ^E'^c CH't C C°'^^°. 

Note that E^ is a Banach algebra for s > [{d+ l)/2]. 

By Proposition 12.121 and Remark 12.141 we immediately obtain: 
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Proposition 2.15 (Tame estimates for second derivatives). Let bo = [{d+l)/2] and suppose s >0 
lies in some finite interval. Then 



(a)|4'([/)(t/^C/^)b. < 



7,T 



(6)(s:([/)({7^f/^)),,,,T< 

c{{u),,,t) {{uis,y,T{uX,T + {u')s,j,Tmbo,T + e{u)s,,,T{uiwMuXa 

In order to obtain a tame estimate for tlie linearized system suitable for Nash-Moser iteration, 
we must recast estimate (j2.57p without the L°° norms of U"' and U^ on the right. First of all, we fix 
the paramater K > 0. For instance, one may take K = 1. This choice is arbitrary because we are 
interested in a small data resullo. We then choose constants So{K), joiK) as in Proposition 12.111 
so that the the estimate (j2.57p holds for s € [0, a], where ce is defined in (|5.60p . For the remainder 
of Section \2. 31 and in Section \5. SX the parameter K is fixed, and 7 is also fixed as 'j = ^q{K). 

Let 

n ■= \U\o,^lo,7,T + {U\xa=o)fMo,'r,T where fiQ ■= [{d+ l)/2] + 2. 

Applying (|2.57p with s = fiQ — 2 we obtain for < e < Eq: 

(2.58) \U\oo,fj.o-2,'y,T + \U\o,fj,o-l,'y,T + {U\xa=o) fio-i,1,T 

< CiK, 7) (|/|o,;.o,7,T + (5)/.o,7,T + {{% + (U)*) ^) ■ 

By Remark 12. 141 if k is chosen small enough, we can absorb the last term on the right in (j2.58p and 
obtain with a new constant C: 

(2.59) |f7|. + (f/U,=o)* < C (|/|o,;.o,7,T + (5)mo,7,t) • 

Substituting (j2.59p in ()2.57p . we have proved 

Proposition 2.16 (Tame estimate for the linearized system). Let K and 7 = ^{K) he as fixed in 
Proposition \2.11\ and suppose \eddU\(jQ,MQ~i + |f7|(^o,Mo < K for e G (0, 1]. Let ^q = [^2^] + 2 and 
s E [0,(5;], where a is defined in (|5.60p . There exist positive constants Ko{'y,T), eq, and C such that 

\U\o,fj.o,y,T + (f/Ud=0>Mo,7,r < Ko, 

then solutions U of the linearized system (|2.54p satisfy for < e < Eq: 

If^b^j, + {U\xa=o)s+l,'y,T 

< C [\f\o,s+2,'Y,T + {g)s+2,'Y,T + (|/|o,Mo7,T + (g) fio ,1 ,t) {\U\o,s+2,'y,T + {U\xa=o)s+2,'y,T)] ■ 



^^If we were interested in a small time result for a given source term G, one would need to fix the constant K in 
terms of G and the parameters 7, T would be chosen accordingly. 
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3 Profile equations 

3.1 The key subsystem in the 3x3 strictly hyperbohc case 

To simplify the exposition we first treat tlie case of a 3 x 3 strictly hyperbolic system and a 
boundary frequency /3 for which there is one single resonance in which two incoming modes interact 
to produce an outgoing mode. This case already contains the main difficulties and is exactly the 
one we emphasize in the example of Appendix [Bj We shall explain later on the relatively small 
changes needed to treat the general case of systems satisfying the assumptions of Section II. 1[ 
The leading profile is decomposed as 

(3.1) V°(x, 61,62, 6s) = ai(x, 61) n + <y^{x, 63) rg 

where (/>2 is the outgoing phase and the resonant triple (?7-i,n2,n3) G Z"^ \ {0} satisfies 

(3.2) ni(t)i=n24>2 + n-i4>s. 
We can thus write 

(3.3) V°„c = o-i(a;,6'i)ri + cr3(2;,6'3)r3, V^^^ = r2(x, 6^2) r2. 
Furthermore we have 

(3.4) V°„cUd=o,0i=e3=eo = (^i^' ^ ^0) e = a{x' , 6q) (ei + 63) , 
so (recall that e = ei + 63, where e^ S spanjrj}, spans ker B n W{f3)) 

(3.5) o-j(x', 0, 6*0) rj = a(x', 6*0)64, i = l,3, 

which determines the trace of Uj in terms of a. 

Applying the operators Ei for i = 1, 3 to (ll.42p (a) and for z = 2 to (ll.42p (b) and using Corollary 
11.261 for ()1.42p (c). we obtain the following system for the unknowns (ui , r2 , (T3 , a) , where A{x',6q) 
denotes the unique function with mean zero such that dg^A = a: 

X^^ai + C10-1 =0 
Xtp^as + csa-i = 

(3.6) '■2^ 
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X,j>2T2 + C0T2 + C2 / (Ti,„i ix, — 62 + — 63] a3{x, 63) d6: 
Jo \ ni m J 

Xlop A + C4^A + C5 ■T2\xa=o + cq (o^)* = -b ■ G* on bi^T , 

where the first three equations hold on i^T, and the constants q are readily computed real constants. 
Here ai^mix, 61) is the image of the function cti under the preparation map 



(3.7) ai{x,6i) = Y,fk{x)e''''^ -^^fkn,ix)e' 



ikniBi 

kei fcez 

a map designed so that the integral in ()3.6p picks out resonances in the product of ai and 173 o 



^Interaction integrals like the one in p.6|) are discussed further in [CGWll] . Proposition 2.13. 
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Differentiating with respect to ^Oi we rewrite the last equation of (j3.6p as 



(3.8) 



Xlopa + C4 a + C5 dggT2\xa=0 + C6 dgg (a 



,G* =: g on bi^T- 



We now set V := {ai{x,9i),a3{x,93),T2{x,62),a{x' ,9o)) and define the interior and boundary- 
operators for the leading profile system: 



(3.9) 



£{V) :- 



X^^ai +ci(Ti 

-^^03 ^3 + C3(T3 



\ 



V 



X^.2T2 + Co T2 + C2 /n " (Ti,„i X, 6*2 H 6*3 ) 0-3(x, 



>d9:, 



B{V) := XLopO. + C4 a + C5 dg^^T2\xa=o + cq (9e„ (' 

In this notation the profile subsystem becomes 

C{V) = in [7t , 
(3.10) B{V) =gm hQx , 

V = Omt<0, 



where the additional relations (j3.5p hold giving the traces of ctijCTs in terms of a. The following 
existence result for the key subsystem is proved in section 15.11 using the tame estimates derived 
below in section [321 

Proposition 3.1. Fix T > 0, let qq := [— 2~] "'"-'-' '^ '■— 2ao + 4, a := 2a — uq, and suppose 
g G H°'^'^ (bQx) ■ Rewrite V as V = {V ,a). If {g)a is small enough, then there exists a solution V 
of the profile subsystem (f3lTO with V G H°'-\VLt), (y'U^=o,a) G H"'-\bnT). 

Remark 3.2. 1) Although the original problem is semilinear with a nonlinear zero-order boundary 
condition, the profile system p.9p has a quasilinear first-order boundary operator and an interior 
operator that includes a nonlinear, nonlocal, integro-pseudodifferential operator given by the inter- 
action integral. The nonlocality arises both from the d^3-integration and from the pseudodifferential 
operator cJi — >• ai^m- 

2) Attempts to solve the system (|3.1Up by a standard Picard iteration lead to a (fatal) loss 
of a derivative from one iterate to the next. The reason is that cJi and a^ have the regularity 
of a (incoming transport equation), and therefore T2 has the same regularity as a. However, the 
equation for a involves the derivative dgf^T2\x^=o and this term induces the loss. Thus, we shall use 
Nash-Moser iteration to prove Proposition 13. 1[ 

3.2 Tame estimates 

With V = (cTi, 0-3, r2, a) and V = (cti, CJ3, f2, d), we compute the first derivatives of £ and B: 



I 



(3.11) 



(a) £'(y)y 



X^^ (Ti -I- c\ a\ 
X^^&s + C3CJ3 



\ 



a;,^^2 + — ^3) <J3{x,93)d93 
ni ni 



X^^f2 +CQf2+ C2 /g "" cri,„i 

+C2 /o^'' cr3,„3 ix, -92 + — 9i] ai{x,9i)d9 

(b) B'(y)V = XLopd + C4d + C5 dggf2\x^=o + 2c6 (a^egd + ddof^a). 
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Here we have used the property 



(3.12) / a3,njx,-—e2 + —ei]&i{x,()i: 
V "3 ni 



dOi = / &i^ni ( X, —02 + —03 ) (Tsix, 6*3) ^6*3, 

Jo V "■! ni J 



which follows readily by looking at the Fourier series of the factors of the integrand. For the second 
derivatives we obtain 



(3.13) 



£"(y)(y^O = c2 



/ 





\ 



+ C ^t, (x, -O2 + -0^ a!(x, 03) d03 

B'\V){V\V'') = 2c6 (d'^Seod' + ii'' de.a'^). 

Proposition 3.3 (Tame estimates for second derivatives), a) Let bo be the smallest integer greater 
than —2— and let s > 0. We have 



C"{V){V^, V%,^ < C \V%,, \VX + \y%n l^lfeo , 



(3.14) 



where C is independent of V , 7, and T. 

b) Let Co be the smallest integer greater than -^ — \- 1 and let s >0. We have 



(3.15) 



(s"(y)(y^ v''))s,, < c (y«).+i,, {v'),, + {v')s+i,, {V^) 



CO ; ' 



where C is independent of V , 7, and T . 

In both estimates (j3.14p . (j3.15p . the constant C can be chosen independent of s in any fixed 
finite interval. 

Proof, a. Moser estimates imply 

(3.16) \al^^(^jc,'^^02 + '^^03)d\{x,03)\Hs(^,^e,) <C {\al^^^^^^^ 

since cig is independent of 02- We have 

The estimate (I3.14p now follows by Sobolev embedding and the fact that 

(3.18) \&lni\s,-y < \^l\s,T 

b. Again Moser estimates imply 

(3.19) {a''de,d')s,^ < C ((d^),,^ {89^)1^ + (d")^- {deod')s,^) , 

so the estimate ()3.15p follows by Sobolev embedding. D 



44 



Next we derive tame energy estimates for the linearized problem 

C'{V)V = f inQT 

(3.20) B'{V)V = gm bflr 

V = Oint<0, 

where / and g vanish in t < 0. We begin with a simple proposition: 

Proposition 3.4. 1. If the phase (pp is incoming, solutions of 

(3.21) ^(jip'^p + (^pf^p = h in Qt, CTp = in t < 
satisfy for 7 large (depending on Cp): 



(3.22) V^Kls,^ < C {ap)s,^ + 



\h\s,^ 
2. If the phase (pp is outgoing, solutions of (|3.2ip satisfy for 7 large (depending on Cp): 



(3.23) \/7kp|s,7 + (o-p)s,7 < C 



\h\s,'y 



V7 
3. Solutions in ujt of 

(3.24) Xioptt + 040 + 2ce{adegd + adg^a) = g, a = in t < 

satisfy for Ck, 7 > 7a' (where K = (a);^i,ooJ; 

C 

(3.25) V^{a)s,-y < —= ((5)s,7 + (a)s+i,7(a)iyi,°°) ■ 

The second term on the right in (I3.25P does not appear in the s = estimate. 

Proof. 1. To prove (I3.23P with s = 0, one considers the problem satisfied by e~'^^ap, multiplies the 
equation by e~'^*ap, integrates dxdOp on Oy, and performs obvious integrations by parts. One then 
applies the L^ estimate to the problem satisfied by tangential derivatives l^~^^^d^i g CTp, |/3| < s. 
Normal derivatives are estimated using the equation and the tangential estimates. The proof of 
(I3.23P is similar. 

2. The proof of ()3.25p is similar, but in the higher derivative estimates one now has forcing 
terms that are commutators involving a. The commutators are linear combinations of terms of the 
form 

(3.26) 7'~""(5f/eo«)(^x%o^^o«) where |/3i| + |/32| = |/3|, |/3i| > 1, 
or linear combinations of terms of the form 

(3.27) ^"^^\d%^^){d^',eode,a) where |/3i| + ^ = |/3|, ^ > 1, 

Applying Moser estimates to p.26p after writing d"^a = d^^^da, we obtain 

.g 28^ (7'""^'(5f'eo")(^<eo^^oa))o,7 < C {{da)Loo{d0^d)m-i,y + {da)rr,-i,j{deod)L--) 

< C {{a)iYi,o^{d)s,-y + (a)s,^(d)t^i,oo) . 
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The factor Ck/ ^/^ on the forcing term in the Lp' estimate allows the first term on the right to be 
absorbed by taking 7 large. 

The estimate of (|3.27|) is similar, but we do not split the d^^'^ derivative, and after absorbing a 
term we are left with -^ {a) l°° {0) s+i,-^ on the right. D 

We now use Proposition 13.41 to estimate solutions of the linearized problem (j3.20p by treating 
the interaction integrals in (j3.1ip (a) and the term c^dQ^T2 in p.lip (b) as additional forcing terms. 
Setting 

(3.29) Vinc,n '■= {(^l,nn CTS.na), Vine = (o"l,Cr3), Vout = T"2, 

estimating interaction integrals as in (I3.16p . (j3.17p . and using (j3.18p . we obtain immediately 

C / 

V Tl '^out\s,'y ~r \ Vout) s,-f _; ^ I 1/ U,7 ~r | yinc,n\L°° \ ''mc|s,7 ~r | i'mc|s,7l Vinc\L° 

(3.30) _ . C ( 

V7l^eKiMtU7 + {deoVout) s,-y < —^ \\def\s,-y + \deVinc,n\L°°\Vinc\s,-y + \deyinc\m,,^\Vinc\L°° ) 

and 

\/uyinc\s,-y S ^ I \Vinc)s,-i H 'J=' 

(3.31) ^ ^^ 



Ck 



V7(^mc)s,7 < -^ \{g)sa + (f^6»oK«t)s,7 + {Vine) s+la{Vinc)\ 



^/7 
This leads to the following "pre-tame" estimate. 

Proposition 3.5. Let //g = [2^] + ^' fi'^ -^1 ^ ^' ^'^^ suppose \Vinc\mi < -f^iO For s > in any 
fixed finite interval, there exist constants C{Ki), 7(i^i) such that for 7 > 7(^1), solutions of the 
linearized problem (|3.2Up satisfy 

\/^\yout,d0Vout,Vinc\s,-y + {Vout , dogVout) s ,-y + ^/^{Vmc)s,■y < 

(3.32) C{Ki 



^— ( l/U+1,7 + \5/s,7 + |^nc|s+l,7l^nc|L°° + {Vine) s+1,^ {Vine) W ''- '°° j 



^/7 
Proof. We add the estimates (I3.30p . (I3.3ip and absorb the terms 

[0.00) -= [\Vinc)s,j 1 xOg^Vout) s,^ ~r \Vine,ni C'gVinc,n\L°°\''ine\s,'Y 

by taking 7 large, after observing that 

(3.34) \Vine,n, dgVinc,n\L°° < C\Vinc,n\fio — C\Vine\fj.o and K = {Vine)w'^-°° ^ C\Vine\no- 



D 



In this Proposition /io ~ [|] + 2 would work, but we make the above choice so as not to have to redefine /io later. 
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and choose constants C{Ki), j{Ki) as in 
^^|. In the remainder of section \3.^ and also 



We now set a := 2a — ao as in Proposition 13.1 
Proposition 13.51 corresponding to the interval s £ [0, d 
in section \U71[ 7 is fixed as j = 'y{Ki). 

To obtain a tame estimate we need to remove the terms depending on Vine on the right side of 
^M- Let 



Applying (|3.32p with s = fiQ — 1 we obtain 
(3.36) 

VtI ^mc I /io- 1,7 + \/7(^mc)^o-l,7 — 1=~ l/Uo,7 + (5)mo-1,7 + ( l^mc|L°° + (Vinc)vi/i.°° j -^2 

By Sobolev embedding if K2 = K2{^,T) is chosen small enough, we can absorb the last term on 
the right in ()3.36p and obtain with a new C 

(3.37) |^mc|L°° + (^mc)Tyi.°= < C (|/|^o,7 + (5)^10-1,7) 

For 7 fixed as above, setting \U\s,'y = \U\s now and substituting (I3.37P in (j3.32p . we obtain the 
estimate in the following Proposition. 

Proposition 3.6 (Tame estimate for the linearized system). Let fiQ = [—5—] + 2 and s € [0, a]. 
There exists k = ^(7, T) > and a constant C depending on k such that if 

(0.00) l^incl/io ~^ \^inc) fiQ < ^) 

then solutions of the linearized system (|3.20p satisfy 

(3.39) \VU + {V), < C [|/|,+i + {g), + (|/U, + (5)^0-1) (I^U+i + (^).+i)] • 

Proof. We have proved the a priori estimate (I3.39P for sufficiently smooth solutions of the linearized 
system. The existence of such solutions now follows by standard arguments, which we summarize 
here for completeness. 

The unknown in the linearized system (j3.20p is ((Xi, (T3, 'r2, d). We can solve the linearized system 
by putting the terms that involve 96)Qf2 or dd^a on the right and replacing the operator dg^, when 
it acts on those terms, by a finite difference operator dg : 

X^^di + cid^ = /i 
X03CJ3 + 030-3 = /2 



(3.40) 



X^,f!i + cof^ = /3 - C2 / " ai,„, (x, ^02 + —93)^lix, e^)de^- 

C2 / \z,nAx,-—02 + —eM{x,e{)ddi 

Jo ^3 m 



Xiopa^ + Cid!" + IcQd^de^a = g- c^d^j!^ - 2cQad^^a!' . 



^*The choice of a is explained in section [5TT1 
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For fixed h G (0, 1] we can solve this system by Picard iteration, where n-th iterates appear on the 
right and (n + l)-st iterates appear on the left. All iterates are in t < and the iterates with 
index zero are all 0. 

We then need an estimate that is uniform in h. This can be done by repeating the existing 
proof of tame estimates, using the operator dg in place of ^g^-^. This gives an estimate like (|3.39p : 

(3.41) iv'^u + {V^), < c [|/i,+i + {g), + (i/u„ + (5)^0-1) {\y\s+i + {y)s+i)] , 

where V^ := {a^ , &i^ , f2 , a^) and C is uniform for h G (0, 1]. Passing to a subsequence, we obtain 
the desired solution of the linearized system. D 

Remark 3.7 (Short time, given data). For a given T > let Ki and 7 = 7(i^i) be as in Proposition 
As we saw above to obtain a tame estimate we need to take |Vmc|^o "I" O^inc)^io small. In our 



formulation of Theorem 11.141 T is fixed ahead of time and we achieve (I3.38P by taking G small in 
an appropriate norm on Qt- For a given G as in (|1.2p vanishing in t < 0, another way to proceed 
is to shrink T; that is, to work on Q^i where < Ti < T is chosen so that 71 := 1/Ti > 7(i^i) and 
so that 

(3.42) \yinc\Hi;0(nT^) + {Vine) Hijo (ujT^) 

is small enough to absorb the terms involving Vine on the right in (I3.36p . One again obtains an 
estimate of the form (|3.39p . where now 

(3.43) \U\s:=\U\H^^^nr,)- 

The iteration scheme described in section 15.11 applies with no essential change to this situation as 
well. 

3.3 The key subsystem in the general case 

Recall that {1, . . . , M} = OUl, where O and I contain the indices corresponding to outgoing and 
incoming phases. We further decompose O = OiL) O2, where Oi consists of indices m such that 
(j)m is part of at least one triple of resonant phases with the property that the other two phases in 
that triple are incoming. For a given m E Oi the phase (pm might belong to more than one such 
triple. 

Now instead of ()3.3p we have 

(3.44) Vine = '^'^(^m,kix,Om)rm,k and Vlut= X] y^^Tm,kix,Om)rm,k, 

meXk=l mSOi fc=l 

since terms T^^k in the expansion of V^^^ vanish if m G O2 as a consequence of (|1.36p and V*^ = V^^^^. 
Recalling that 

(3.45) ^ = X^ X^em,fc, where em,k e span{r„,A:}, 

mGX fc=l 
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we see that in place of (|3.4p we now have 

(3.46) V°„cUd=0; 9m=eo,rn(ii = a{x' , 6*0)6 = ^ ^ aem,k = X] X] '^rn,k{x' , 0, Oo)rm,k, 

meX fc=l mel A;=l 

and thus 

(3.47) am,k{x',O,0o)rm,k = a{x',0o)em,k iov m£l, k = l,...,Uk^. 

Next we derive the formulas for C{V) and B{V) in the general case. The unknown is now 

(3.48) V = {ajn,k,meI,k = l,...,Uk^; Tra,k,m e Oi,k = 1, . . . ,Uk^; a) 
Suppose q G Oi and that {(t)p,4>q,(l)s) is a resonant triple such that 

(3.49) np(j)p = UgCpq + ns(j)s where p,s £ I and gcd(np, Uq, Ug) = 1. 

Applying the projectors Em,k, m £ I, k = 1, . . . , v^m to (11.42j) (a) and the projectors Eq^i, q £ Oi, 
I = 1,... ,uu to (ll.42p (b) we obtain 



CiV) 

(3.50) I ^ ^^I'fep sr'^ks jk,k' r2TT 

1 2^A:'=1 "'q,l Jo 

(similar ); Q' G Oi,/ = l,...,i/fc^ y 



^4>m^m,k + Cm,k(^m,k', m £l,k — 1, . . . , V^^ \ 



Here "(similar)" denotes a finite sum of families of integrals similar to the family given explicitly 
in (|3.50p . Here "family" refers to the sum Ylkk'- One such family corresponds to each distinct 
resonant triple involving the outgoing phase (pq and two incoming phaseso The values of the real 
constants Cm,k, d \ are not important for our analysis, but for example the d '^ are given bjCj 

(3.51) d^^'f = —iq^i ■ [d^D{0){rp^k,rs,k') + dvD{0){rs,k',rp^k)]- 

By a computation similar to the one that produced ()3.8p we obtain from ()1.42p (c) 

(3.52) B{V) = XLopa + fia+^Yl fi^i^e.Tq^i + f2de,{a^). 

qed 1=1 

for some real constants /i, /2, fq^i- For example, we have f2 = —b • [tp'{0){e,e)]. Thus, the system 
()1.42p may be rewritten 

CiV) = in r^T 

(3.53) B{V) = -b • de^G* := g on b^r 

y = in t < 0, 



29 



We do not distinguish between [cj>p,cj>q,(j)s) and [(j>p,(j>s,(j>q). We do distinguish between {(j)p,(j)q,(j}s) and 

i>p,(j>q,cl>t)- 

^"We have suppressed indices r, s on the d \ . 
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where the relations ()3.47p hold. 

It is now a simple matter to write out the expressions for the first and second derivatives of 
C and B. For example, just as the interaction integral in (|3.9p gave rise to two integrals in the 
expression ()3.1ip for CJ (y^ in the 3x3 case, it is clear that each integral in ()3.50p will give rise 
to two integrals in the new expression for CJ{y\ The tame estimates for second derivatives are 
proved exactly as before, and Proposition 13.31 holds verbatim in the general case. Proposition 13.41 is 
used exactly as before to prove estimates for the linearized system. With the unknown V as given 
in p.48p and after defining Vine-, Vout, Vine, Vout in the obvious way, we see that the "pre-tame" 
estimate of Proposition 13.51 and the tame estimate of Proposition 13.61 hold verbatim in the general 
case. The iteration scheme of section [5T] depends only on the tame estimates. Thus, it applies here 
without change and Proposition 13.11 holds verbatim in the general case. 

Once the key subsystem is solved, we can easily complete the solution of the full profile system 
()1.35p . ()1.36p . The precise result for the full system is proved below in Theorem 15.111 
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4 Error analysis 

Here we carry out the error analysis sketched in section 11.51 beginning with the proof of Propo- 
sition [T^Sl 

Proof of Proposition \1.2R 1. Noncharacteristic modes. We write 

M 

F(x, 9) = Fo{x) + Y, Fa{x) e^"-' + Y. Y. ^"(^) ^"■'' 

a^e rn=l aeCm\{0} 

and recall that the sums are finite. Set 

M 
na = Y "i ^^^^ ^ = (^li • • • i^Af)- 

Since EF = 0, we first note that Fq vanishes. For any a, there holds 

and when a ^ C we look for Ua{x) such that 

(4.1) £o(9eo, %)f/a(2:)e^'^"^o+*("-^)«'' = F„(x)e*""^«+'^("-^)«^. 

This holds if and only if 

iL{na/3,a ■ u}_)Ua = Fa- 

The matrix on the left is invertible so we obtain a solution of (|4.1|) for a ^ C. 
2. Characteristic modes. When a E Cm \ {0} we have a • w = ?^o w^, so 

We can write 

aec„A{o} A;eZ\{0} 

where 

{aeCm\0,na=A;} 

Since -Em-^ = 0, we have for each A; G Z \ {0} that PmFm,k{x) = 0, so now we look for Um,k{x) 
such that 

Co{de,,d^,)Um,k{^) e^'^(^«+^'"«^) = (/ - Pm) Fm,k e'^^'^^+^-^^l 

The latter relation holds if and only if 

iL{k/3, ku^) Um,kix) = ikL{d4)m) Um,k{x) = {I - Pm)Fm,k{x), 
which is solvable even though L{d(j)m) is singular. Finally, we take 

M 

U{x, 00, ^d) = Y ^-(^) e^"'^^"+^("-^)«^ + E E ^-.'^(^) 6^'^'°+-'"^''^ , 

ai^C in=l fceZ\{0} 

which solves (|1.46p as claimed. D 

51 



The existence theorems for profiles and for the exact solution to the singular system, Theorems 
15.111 and 15.131 respectively, are stated and proved in section \5\ we shall only use the statement of 
these theorems here. In order to formulate the main result of this section we must make some 
preliminary choices. 

Choice of a and a. The conditions on the boundary datum G{x' ,6q) are slightly different 
in Theorems 15.111 and 15.131 We need to choose a, a, and G{x',6o) so that both Theorems apply 
simultaneously. We also need a large enough so that we can apply Proposition 12.91 in the step 
(j4.24p of the error analysis below. These conditions are met if we take 

(4.2) a = max (d + 9, [{d + l)/2] + Mq + 3) and a = 2a-[{d + l)/2] 



and choose G E H°'{bVLT) such that (G) jja+2 (^f,^^^ is small enough. As in Theorem 11.141 we now 
rename the numbers in (j4.2p as a and o. Applying Theorems 15.111 and 15.131 we now have for 
< e < £0 an exact solution Ue{x,6Q) G E°'~^{Q.t) to the singular system (ll.lSp and profiles 
V^{x,e) G H^-^inr), V^{x,e) G H^-'^inr) satisfying the equations (fL35]) and (fL36]l . 

Approximation. Fix 6 > 0. Using the Fourier series of V^ and V^, we choose trigonometric 
polynomials Vp{x,6) and Vp{x,6) such that 

(4.3) |V° - VpV.-i(n^) < 5, \V' - V^V-^cn^) < S. 

We can smooth the coefficients so that Vp and Vp lie in H'^{Qt) and so that (j4.3p still holds. 
Having made these choices, we can now state the main result of this section, which yields the final 
convergence result of Theorem 11.141 as an immediate corollary. 

Theorem 4.1. We make the same Assumptions as in Theorem, \1.14\ ^nd let a and a be as just 
chosen. Consider the leading order approximate solution to the singular semilinear system (jl.lSp 
given by 

(4.4) Z^°(x,0o):=V°(x,0)|,_(,^^^), 

and let Us{x,9o) G E^~^{Q,t) be the exact solution to (jl.lSp just obtained. Then 

(4.5) liui\Ue{x,eo) -l(^{x,eo)\Ea-^nr) = 0- 

6—>-0 ^ ' 

The following lemma, which is proved in [CGWll] . Lemmas 2.7 and 2.25, by a simple argument 
based on Fourier series, is an important tool in the proof. 

Lemma 4.2 (Relation between norms). For m G N suppose f{x,9j) G H^'^^(Q,t), ind set 
fe{x,eo) = f{x,eo + Ujf). Then 

(4.6) \fe\E^<C\f\Hm+i^^^y 



Proof of Theorem \4-l\ We shall fill in the sketch provided in section 11.51 
1. First we use Proposition 11.281 to construct h(p{x,9o,^d) satisfying 

(4.7) Co{de,,d^,)U^=[-{I-E){L{d)V'p+DiO)Vl + d,DmvlV^p))] \e^^eo,i. 
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The function U"^ is a trigonometric polynomial of the form (jl.47p with H°° coefficients. We then 
define the corrected approximate solution 



(4.8) 



Ue{x,ea) := {v\x,e) + ev\x,e))\,_^^,^^^_,) + e'ulix^e^,'^). 



Since V^ £ H''-'^{^t), Lemma SJ implies Ue € ^''-^(Ot). 

2. Next we explain (I1.48P and make precise the norms used on the right there. Using the 
identity ()1.44p we compute 



01,o^ 



(4.9) h,{Ue) = e (>Co(9eo,%)^p)lg,=£^ + {L{d)V^ + D{Q)V^ + d,D{^)V''V 



le^{0o,^) 



+ 0{e^). 



Here the profile equations ()1.20p (a).(b) imply that the terms of order e ^ and e^ vanish. Using 
(j4.7p we can rewrite the coefficient of e in (|4.9p as 



(4.10) [L{d){V' - VD + Dmv' - VD + a,Z)(0)(V°V° - V°V°)] |,^(,^,_^) 

+ [E {L{d)Vl + D{^)K + 5.^(0)V°V°)] le^(,„,f.) 



A + B. 



Using (j4.3p . Lemma W?2\ and the fact that E'^{Qt) is a Banach algebra for s > [{d + l)/2], we see 
that 

(4.11) |^|£;-4(o^) < K6. 
To estimate B let 

(4.12) F = L{d)V^ + D{0)V^ + 9^L>(0)V°V° and Fp = L(9)VjJ + D{0)Vl + 9^,F'(0)V°V°. 



The profile equation (ll.SGp fbl implies EF = 0. Using continuity of the multiplication map (|1.25p . 
we see that (j4.3p impliecH 



\F - Fp\j^a-3;2 < K5. 



(4.13) 

From the continuity of E : H^ — t- H^ and Lemma 14.21 we then obtain 



(4.14) 



\B\E^-*{nT) 



{EF,) 



5-^(^0,^) 



£"-4(0^) 



(F(F-Fp))|,^(,^,^,, 



E'^-^(nT) 



<K6. 



3. The O(e^) terms in ([0|) consist of 



(4.15) 



£2 (L(a)Z^2^x,^o,Cd)) 



9o,^) 



£"-4(f7r) 



<e'C(5), 



as well as terms coming from the Taylor expansion of D(eZ//e)ZYe) like {e^dvD{0)V^V^] 
of which satisfy an estimate like (j4.15p . Setting Ri;{x,9q) := I^si'^e)^ we have shown 



all 



(4.16) 



\Rs\E'^-Hnr)<<K^ + CiS)e). 



Here H!^ ' denotes the space defined in (|1.24[) . but with the obvious restriction on the domain of t. 
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4. The boundary profile equations (jl.22p and the fact that the traces of V'' and V^ He in 
H"-^^{bQ.T) and H°-~'^ {bQ.T) , respectively, imply 

(4.17) (r-e(x',eo)>^.-2(,f,^) < C{5)e\ where r, := i;{eU,)U, - eG{x' ,60). 
Indeed, these O(e^) terms include 

(4.18) {e^BU^pix',0,eo,0))j^^,,^^^^^ < C{5)e\ 

and other terms satisfying the same estimate coming from the Taylor expansion of ip{ehlf.)'Ue- 

5. Next we consider the singular problem satisfied by the difference W^ := U^ —Ue- 



(4.19) 

where 
(4.20) 



ddWe + A (d^, + ^^ We + D2{eUe,eUs)We = -Re 

ii2{eUe,eUe)We = -r^ on Xrf = 
We = in t < 0, 



D2{eUe,eUe)We := D{eUe)Ue - D{eUe)Ue = 

D{eUe)We+(f d,DieUe + se{Ue-Ue))ds] {Ws,eUe) 



and 'il)2{£U£,£hle)We is defined similarly. Since U^ G E"" ^(f^r) and 14^ € E'^ ^{^t) a short compu- 
tation shows 

(4.21) 

^lj2{eUe,eU,)We = tl^(eUe)W, + d,i^{eU){W„eU,) + 0{C{6)e^) = B{eU,eU)W, + 0{C{5)e^), 

where the error term is measured in H°'~^{hQ.T) and B is defined in (II. 9p . Similarly, 

(4.22) D2{eUe,eUe)We = V{eU,eU)W, + 0{C{6)e^) in ^"-^(fiT). 

Thus, using ()i36]) and KT7\i we find 



(4.23) 



ddWe + A ( 5,, + ^^ We + V{eUe,eUe)We = e{K6 + C{6)e) in E'^-^rtr) 

-2\ :„ T_ra— 3/ 



Bi£Ue,eUe)WeU,=o = 0(C(5)e^) in H^-'\b^T) 
We = in t < 0. 

Applying the estimate of Proposition 12.91 we obtain 

(4.24) \We\E^^^^)<K5 + C{5)e, 
which implies 

(4.25) \U, - ^^°bo(f,^) <K5 + C{5)e. 
Fixing first 6 small and then letting e — )• we have shown 

(4.26) lim|[/e-W,%o(n^)=0. 

The family Ue —Us, < e < eoi is bounded in E°-''^{Q.t), so by interpolation (|4.26p implies 

(4.27) lim|^e-^^°b-3(o^)=0 

as required. D 
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5 Nash-Moser schemes 

5.1 Iteration scheme for profiles 

A good reference for the Nash-Moser scheme is |AG07j . The method depends on having a 
family of smoothing operators with the following properties. For T > 0, s > 0, and 7 > 1, we let: 

(5.1) F^'^i^r) ■■= {u G H^{nT),u = for t < 0}. 

Lemma 5.1 ( |Ali89| . section 4). There exists a family of operators Se : F!^{^t) — ^ n^>oi7 (^t) 
such that: 

(a) \Seu\i3 < Ce^'^-''^+\u\a for all a, (3 > 

(5.2) {b)\Seu-uy<C6^^-'''>\u\a, < /3 < a 

(c) l^zSeulp < Ce^^-'^-^^\u\a for all a,/3 > 0. 

The constants are uniform for a, 13 in a bounded interval. 

There is another family of operators S0 acting on functions defined on the boundary and satis- 
fying the above properties with respect to the norms {u)s on bfl 
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Tl 

Description of the scheme. Our goal is to solve the problem (|3.10p : 

C{V) = in ^T 

(5.3) B{V) = g in bQr 

V = Oint<0. 

The scheme starts with Vq = 0. Assume that I4 fire already given for k = 1, . . . ,n and satisfy 
Vfc = for t < 0. We define 

(5.4) Vn+l = Vn + Vn, 

where the increment Vn is specified below. Given 9q> 1, yjq set 9^ '■= {Oq + n)^'^ and work with 
the smoothing operators Sg^. We decompose 

(5.5) C{Vn+l) - C{Vn) = C'{Vn)Vn + c'^ = C {SeJ^nWn + e'„ + e'^, 

where e^ denotes the usual "quadratic error" of Newton's scheme and e^ the "substitution error" . 
Similarly, 

^((K+i)u,=o) - ^((K)u,=o) = ^'((K^)u,=o))(KU,=o)) + e; = 

(5.6) 

B' {{S0^Vn)\xi=o){Vn\x^=o) + &„ + e". 
The increment Vn is computed by solving the linearized problem 

C'{Se^Vn)Vn = fn 

(5.7) B'((5,„K)U,=o)(KU,=o) = gn 

T4 = in t < 0, 



^For u defined on Q,t we do not necessarily have equality of {S0u)\x^=o and Se{u\xj^=o) 
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where /„ and gn are computed as we now describe. 
We set e„ := e^ + e'^ and e^ := e'„ + e". Given 

yo:=0, /o:=0, go := Se,g, Eq := 0, Eo := 0, 

(5.8) Vi,...,Vn 

fi,- ■ ■ , fn-i, gi, ■ ■ ■ ,gn-i, eo,...,e„_i, eo,...,en_i, 

we first compute for n > 1 tlie accumulated errors 

n— 1 n— 1 

(5.9) En:=^ek, En:=^ek. 

fc=o fc=o 

We tlien compute /„ and gn from the equations 

n n 

(5.10) 5^/fc + 5e„i?„ = 0, 5]5fc + ^e„^n = 5e„<7, 

fc=0 fc=0 

solve (|5.7p for y„, and finally compute _Vn_|_i from (j5.4p . 
Next e„, e„ can be computed from ri 

,, ,,, £(K+i)-/:(K) = /„ + e„ 

(5.11) 

B{{Vn+l)\xi=o) - S((K)Ud=0) = 5n + en- 

Thus the order of construction is 

(5.12) 5* (e„_i,e„_i) -^ {En, En) -^ {fn,gn) "^ K "^ K+1 "^ (en,e„) -;>... 

Adding (j5.1ip from to n and using (|5.10p gives 

,,,,, C{Vn+i) = {I-SejEn + en 

(5.13) _ _ _ 
m(K+i)|x,=o) - 5 = (Se„ - /)5 + (/ - 5eji5;„ + e„. 

Since Se,^ — >■ / and 5*51,^ — ;■ / as n — ;■ oo and we expect {en, en) — ?■ 0, we formally obtain a solution 
of (j5.3p in the limit as n — )• c«. 

Induction assumption. Let A^ := 9n+i — On and observe that 

(5.14) ^ < A, = VOl + l- On < 7T^ for all n G N. 

6Un ^(7n 

With //Q = [^2^] + 2 as in Proposition 13.61 we now set ao := //q — 1 and fix a choice of integers 
ao < a < a, whose values are explained below: 

(5.15) a = 2ao + 4 and a = 2a — qq. 
Given 5 > our induction assumption is: 

(Hn-i) For ah fc = 0, . . . , n - 1 and for all s € [0, a] n N 

(5.16) \Vk\s + {Vk)s < SOf-^-^Ak. 



'in the estimates of e„ and Cn, we instead use the formulas (|5.2ip . H5.25|) 
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The main step in the proof of Theorem 15.111 is to show that for correctly chosen parameters 
6 > (small) and ^o ^ 1 (large) and for small enough g, (Hn_i) implies (H„). At the end we will 
verify that (Hq) holds for for small enough g. 

First we state some easy consequences of (H„_i). 

Lemma 5.2. If 6q is large enough, then for k = 0, . . . ,n and all integers s G [0, d] we have 



3(s-a) + 

(5.17) \Vk\s + {Vk)s< ''^" ^ 

I u d log fJk, a = s 



CS0)!-''>+, a^s 



Proof. This follows by writing Vk = Vq + X],=o ^j ^"^^ using the triangle inequality and an elemen- 
tary comparison between Riemann sums and integrals. D 

Lemma 5.3. // ^o 'is large enough, then for k = 0, . . . ,n and all integers s G [0, q + 2] we have 

(5.18) |5.>^.|,4«f7-^» . 

[Cdlog9k, a = s 

Moreover, for k = 0, . . . ,n and all integers s G [0, d] we have 



(5-19) \il - Se.nU < .^^,_^ 



0691-'' log Ok, ^^« 
0591'°, s>a 



Proof. This follows from Lemma 15.21 and the properties of the 3$. For example, we have 

|(/ - Se,)Vk\s < 2|yfc|, < CJ0^-" for s > Q 
^ ' ' \{I- Se,)VkU < Cr-^lFfcU < Cd9'-"log9k for s < a. 

Estimate of the quadratic errors. From (j5.5p and ()5.6p we have 

(a)e'fc = C{Vk+i) - £{Vk) - £'{Vk)Vk = / (1 - r)£"{Vk + rT4)(Vfc, 14)dr 

(5.21) '^\ 

(6)4 = B{Vk+i) - B{Vk) - B'{Vk)Vk = / (1 - r)B"{Vk + TV^m, Vk)dT 

Jo 

where the arguments in (|5.2ip (b) are evaluated at Xd = 0. 

Lemma 5.4. 1) For large enough 9q we have for all k = 0, . . . ,n — 1 and all integer s G [0, a] 

(5.22) \e',\s < C6'9i'^'^~'Ak, 

where Li{s) = s + oq — 2a — 2. 

2) For large enough 9q we have for all k = 0, . . . ,n — 1 and all integer s G [0, d — 1] 

(5.23) {e',)s < C6'9i''^'^-'Ak, 
where L2{s) = s + ao — 2a — 1. 
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Proof. Using (j5.2ip (a). Proposition 13.31 and the fact that ao > Bq, we have 

(5.24) WkU < C\Vk\s\Vk\ao- 

The estimate (|5.22p then follows by applying the assumption (j5.16p and using A^ ~ ^. The 
estimate (j5.23p is proved the same way; the restriction s £ [0, a — 1] reflects the loss of one 
derivative in ()3.15p . 

D 

Estimate of the substitution errors. From (|5.5p and ()5.6p we have 

(a) el = [ C"{Se,Vk + r(T4 - Se.Vk)) (Vk, {I - SejVk) dr 

(5.25) ^° 

(6) 4 = J B"{Se,Vk + T{Vk - Se.Vk)) (t4, (I - Se,)Vk) dr 

where in (15.25p (b) we have, for example, written Sof^Vk for {Sg^Vk)\xa=o- 

Lemma 5.5. 1) For large enough 0q we have for all k = 0, . . . ,n — 1 and all integer s G [0, a] 

(5.26) \el\s < C6'e^,'^'^~'Ak, 

where ^3(5) = s + uq — 2a + 1. 

2) For large enough Oq we have for all k = 0, . . . ,n — 1 and all integer s £ [0, d — 2] 

(5.27) {e'Ds < C6^et'^'^-'Ak, 
where Li{s) = s + oq — 2a + 3. 

Proof. Using (j5.25p (a) and Proposition 13.31 we obtain 

(5-28) le'fc'l, < C (\Vk\s\{I - Se,)Vk\ao + \{I- Se.WklMl 



The estimate (j5.26p now follows from (H„_i) and Lemma l5.3i The estimate ()5.27p is proved the 
same way, after using the trace estimate 

(5.29) ((/ - Se,)Vk)s+i <C\iI- %)T4|.+2. 

The restriction s € [0, d — 2] reflects the subscript s + 2 in (I5.29p . D 

Estimate of {En, En) and {fn,9n)- Since e^ = e'f^ + e'^, and e^ = e'j^ + e'^,, we have 
Lemma 5.6. There exists 9q sufficiently large so that 

(5.30) \En\a < C5^e^'^^^ and {En)a~2 < Cd^et^^-^\ 



Proof. Viewing En = X]fc=o ^*: ^^ ^ Riemann sum and using ^3(0) > OJfj we obtain the estimate 
of En from ()5.22p and (j5.26p . Since L4(d — 2) > 0, the estimate of En is similar. D 



*This determines a in (|5.15ll . 
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From (jS.lOp we have 
(5.31) _ _ ^ -. -. _ _ 

Lemma 5.7. There exists 6q sufficiently large so that for s G [0,(5 + 1] we have 

(b) {gn)s < CSH^^'^'^^n + Cei-^-'{g)^An. 
Proof Using (IS2D(c), dEMJ, and s - d + ^3(0;) = ^3(5), we find 

(5.33) |(5,„ - Se„_,)Er,-i\s <C T e'-'''^\En-iUde < C6^9^t^f-'Ar^. 

From (J5.22p . (j5.26p and the properties of So we readily obtain 

(5.34) |5e„e„_i|, <C<520^«W-iA„, 

and this gives (j5.32p (a). 

The first term on the right in (j5.32p (b) arises similarly. With 

(5.35) {{Se,, - Se^_,)g)s <C T e-'-"-\g)^de < C0r""'(5)«A„ 

we obtain (I02|) (b). D 

Induction step. We claim that for 6 > sufficiently small, the estimate for the linearized 
system (|3.39p applies to (|5.7p and gives for s G [0, a]: 

(5.36) \Vn\s + {Vn)s < C [\fn\s+l + (S'n)s + (|/n|ao+l + {9n)ao) (l-Se^KU+l + (•S'e„K)s+l)] 

Indeed, (|5.18p and a > ao+2 imply that for (5 > small enough, the requirement (I3.38P holds r^l For 
the terms involving /„ and gn, except {gn)ao-, we substitute directly into (j5.36p the corresponding 
estimates from Lemma l5.7i For {gn)ao we have 

(5.37) {gn)ao < C (<52e^(-«)-iA„ + 0-''-\g)^,+^+iAn) , 

where the last term arises from (j5.35p with s = ao and a replaced by a + ao + 1 ■ We also use 

(5.38) {Se,yn)s+i < \Se„Vn\s+2 < 069^:-^^-"^++' 

and a similar estimate for jS'e^Vnls+i, which follow directly from (I5.18p . 
Since ^4(5) > L^i^s + 1) this gives for s G [0, q]: 

(5.39) 

\Vn\s + {Vn)s< 



c 



S2qU{s)-i^^ + 0r"-i((7)„A„, + U2^^("o)-iA„ + 0--2(5)„„+,+iA J 59^^ 



(s+2-a)++l 



^We use a trace estimate like (|5.38p here as well. 
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For s G [0, a] the parameters qq and a (recall (|5.15p ) satisfy: 

L^is) < s — a 

(5.40) L4{ao) + {s + 2-a)+ + l< s-a 

(s + 2 - a)+ < s. 

Thus, we have proved (H„), which is the content of the following lemma. 

Lemma 5.8 (H„). If 5 > and {g)a/^ cLfe sufficiently small and 9q sufficiently large, we have 

(5.41) \Vn\s + {Vn)s < (56l^~""^A„ for all integer s £ [0,d]. 
Still assuming (H„_i) we now show: 

Lemma 5.9. Suppose n>l. If 5 > is sufficiently small and Oq sufficiently large, we have 

(a) \C{Vn)\s < 50'"""^ for all integer s G [0, d] 

(5.42) , 

(6) {B{Vn) - g)s < 5^r"" for all integer s E [0, d - 2]. 

Proof. From ()5.13p we have 

(a) l^(K)|s < |(/ - Se„_,)En-i\s + Wn-ils 

{b) {B{Vn) - g)s < {CSe^., - I)g)s + {{I - Se„_,)K-i)s + (e„-i).. 

Using ()5.2p and the above estimates of Sn-i and e^-i, we find 

\en-l\s < C6X''-'^~'^n, 

which imply ()5.42p (a) since ao + 2 — a < and -^3(5) < s — a. 

The last two terms on the right in (j5.43p (b) are estimated similarly. To finish we use 

(5.45) ((5,„_, - I)g)s < Cet^''-^\g)a-2 for s < d - 2 

and observe that s — d + 2<s — q — 1. D 

We now fix 5 and ^o as above and check (Hq). 
Lemma 5.10. // {g)a is small enough, then (Hq) holds. 
Proof. Applying the estimate for the linearized system to 

£'(0)^0 = 

(5.46) , ^ . 

B'{0)Vo = Se,g 

we obtain for integer s G [0, d]: 

(5.47) \Vo\s + (Fo). < C{Se,g)s < C /^o""^^^- '^" . 

[{g}a, s < a 

Thus, (Ho) holds if {g)a is small enough. 

D 
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We can now complete the proof of Proposition 13.11 
Proof of Proposition l^'.Jl We have 



n— 1 n— 1 



(5.48) Vn = Vn-l + Vn-l = Vo + ^Vk = Y,Vk. 

k=0 k=0 

Let /U = a — 1. Since 6f^ ~ \/k we have by (H„) 



(5.49) ^|T4|^ + ^(T4)^<5^0,'Afe<C^A;-2 <oo. 

jk=0 fc=0 k k 

Thus, for some V as described in Proposition 13.11 Vfe — >• ^ in H'^{Qt) and Vk\xa=o ~^ ^Utj=o in 
H^'{nT)■ This imphes 

(5.50) aVk) -^ C{V) in H''-\nT) and ^(14U,=o) ^ ^(^U,=o) m F'^-1(61]t). 
Applying Lemma fS.QI with s = fi — lwe conclude that F is a solution of the profile system (IS.lOp . D 

Having solved the key subsystem we can now easily complete the solution of the full profile 
system (jl.35p . (|1.36p and obtain the following result. 

Theorem 5.11. Fix T > 0, let uq = [— 2~] + 1, a = 2ao + 4, a = 2a — ao, and suppose 
G G H'^^^{Q,t)- If {G)a+i is small enough, there exist solutions 

(5.51) v° = vL e h^-H^t), v'=y} + vL + vL g h'^-^^t) 

of the full profile system (jl.35p . (jl.36p satisfvinr^l 

(5.52) VL = EVl^ G //"-^(f^T), (i?VDU,=o,e,=eo G H"-\bnT) 

y} G i^"-2(i7T), (/ - E)vi^ G i/°-2(f^T), ^vL G i/°-^(i^T). 

These statements remain true if a is increased and if a>2a — a^. 

Proof After the subsystem pri2|) is solved we know V° = V°„^ = -EV",^, V^ut = EV^ut, and 
a, and these functions have the regularity described in Proposition 13.11 Taking the mean of 
equations ()1.36p (b)(c)(d). using the fact that the mean of the quadratic term in (jl.36p (b) lies in 
H°'~^{Q,t), and applying the result of [CouOSj to the resulting weakly stable system, we conclude 
yi ^ H'^-'^inr)- From (fL36]l fa) we find 

(5.53) (/ - E)V' = {I- E)Vl, G F"-2(f7r). 

It remains to determine EV}^^. Since the solvability condition (ll.41|) holds, we can make a choice 
of EVl^^\xj^=ofij=6o G -ff"~^(6r2T) satisfying the boundary equation (|1.40p . whose right side is now 
known and lies in ff "~^ (br^T-) r^Finally. we determine the components of EVinc by solving the 
transport equations determined by (jl.36p (b). the choice of initial data, and the initial condition 
(jl.36p (d). Observe that the interaction integrals corresponding to the quadratic term in (|1.36p (b) 
lie in H'^~^{Q.t)- □ 



^^Here when we write Vf„;, G H"^^{Q,t), for example, we mean that the individual components of Vf„c lie in that 
space. 

^^All terms on the right in p.4Up lie in H°'^^{hQ,T), except the term involving L{d). That term is actually more 
regular than H"'^ {bOr) , but we do not wish to introduce more refined spaces to capture this. 
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5.2 Iteration scheme for the exact solution 

The Nash-Moser scheme for the exact solution wih use the scale of spaces E^ j, on VLt and 
H^ rp on hVtT- Since T was fixed at the start and 7 was fixed in section 12.31 we now drop these 
subscripts in the notation for norms and function spaces. For s > we let 

(5.54) F'^ := {u{x, Oq) £ E',u = Oiort< 0}. 
Moreover, we shall now denote E'^ norms simply by \U\s and H^ norms by {U)s- 
Lemma 5.12. There exists a family of operators 6*0 : F" —t- n^>oF^ such that: 

(a) \Seu\f3 < C0('^~")+|u|a for all a,/3 > 

(5.55) {b)\Sgu-uy<Ce^^-''^\u\^, 0<l3<a 

(c) \4eS9u\p < Ce^^-'^-^^\u\a, for all a,/3>0. 
du 

The constants are uniform for a, (3 in a hounded interval. 

There is a family of operators Sg acting on functions defined on the boundary and satisfying the 
above properties with respect to the norms {u)s on bil.T, o,nd we have 



(5.56) 



{Seu)[j,^=Q = Se{u\x^=o). 



Proof. Let Sg be a standard family of smoothing operators, for example as in |Ali89] . acting in the 
{x\Oq) variables on the scale of spaces H^ . For U £ E^ simply treat Xd as a parameter and define 



(5.57) 



Sell = SeU{-,Xd,-)- 



The properties (j5.55p then follow immediately from the corresponding properties of the operators 
Se. □ 

To avoid excessive repetition we use the notation and arguments of section 15.11 as much as 
possible, and just point out where changes are needed. Thus, we now denote the solution to the 
semilinear singular problem ()1.18p by V instead of U, write g instead of G, and rewrite (jl.lSp as 



(5.58) 
where 
(5.59) 

We now lelH 
(5.60) ao -- 



C{V) 
B{V) 



C{V) = on ^T 
B{V) =gon bQr 
V = in r < 0, 






V + D{eV)V 



d+1 



ai 



d + 1 



+ Mq, a = max(2ao + 3, ai + 1), a = 2a — oq- 



The main result of this section is the following proposition. 



•^^The parameter a is determined so that L2{a) > for L2{s) as in Lemma 15.181 The definition of a is chosen so 
that ai < a and the conditions (|5.80p hold. 
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Theorem 5.13. Fix T > 0, define a, oq, and a as in (|5.60p . and suppose g G H"^, where g is the 
same as the function denoted G in (jl.lSp . There exists £q > Q such that if {g)a+2 is small enough, 
there exists a solution V of the system (|5.58|) on $7^^ for < e < Eq with V G E°'~^ , V\x^=o S H°'. 
Thus, Us = V is a solution of the singular system (jl.lSp on Qt for < e < eg- These statements 
remain true if a is increased and if a > 2a — a^. 

The hnearized singular problem (I2.54p is now written 

C'{V)V = f onQT 

(5.61) ]3'{V)V = gon bVLr 

y = in t < 0. 

With this notation the description of the scheme in section 15.11 starting at line (j5.3p applies here 
word for word down to line (I5.14p . 

Remark 5.14. (a) In order to apply the tame estimate of Proposition 12. 16] to the linearized system 
(|5.6ip . by Sobolev embedding (Remark 12. 14p it suffices to have 

(5.62) \eddV\a,-i + \VU, < K' for e £ (0, 1], and I^U^+a < k 

for some constant K' depending on K and k. as in Proposition 12.161 In fact we use the slightly 
weaker (because we use E^ norms on the right) estimate for s G [0,(5]: 

(5.63) \V\s + {V)s+i < C O/U+i + {g)s+2 + (|/Uo+i + (5)^0+2) i\U\s+i + {U)s+2)] • 



(b) By Proposition 12.151 when \V\ao < K' , the tame estimates for second derivatives now take 
the form 

{a)\C''{V){V'',V%<c(\V%\vX, + \V%\V''U+e\V\s\V''UvXo 

(5.64) ; 

{b){B"{V){V^,V'))s < C {{V'')s{V')o., + {V')s{V'^)o.,+e{V),{V'^)o.,{V') 



With a and a redefined as in (15.601) . for a given 5 > the induction hypothesis {Hn-ij is now 

(H„_i) For all A; = 0, . . . , n - 1 and for all s G [0, d] n N 

(5.65) \Vk\s + {Vk)s+i < <50r"-'Afc. 

Lemmas 15.21 and 15.31 are now replaced, with no real change in the proofs, by the following two 
lemmas. 

Lemma 5.15. // 9q is large enough, then for k = 0, . . . ,n and all integers s G [0, d] we have 



^(s-a) 

(5.66) \Vk\s + {Vk)s+i< ;■,■" ^ 

[Co log Uk, a = s 



C6el'-''>+, a^s 
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Lemma 5.16. // 9q is large enough, then for k = 0, . . . ,n and all integers s £ [0, a + 2] we have 

C6log6k, a = s 



(5.67) \Se,Vk\s + {Se,Vk)s+i< 



For k = 0, . . . ,n and all integers s € [0, d] we have 



(5,68) i(/ - se,n\. + ((/ - .%,.nu, < ':'%°J°'"'- ' ^ " , 

\CdOl , s > a 

We have used (|5.57p for the estimate on traces in Lemma 15.161 In place of Lemma 15.41 we now 
have: 

Lemma 5.17. 1) For large enough Oq and small enough 5 we have for all k = 0, . . . ,n — 1 and all 
integer s G [0, a] 

(5.69) \e',\s < C6^ei'^'^-'Ak, 

where Li{s) = max(s + ao — 2a — 2, (s — a)+ + 2ao — 2a — 1). 

2) For large enough 9q and small enough 6 we have for all k = 0, . . . ,n — 1 and all integer 

s £ [0, a] 

(5.70) {e',)s+i < C752^^i(^)-iAfe. 

Proof. Again we use the formulas in ()5.2ip . By Lemma 15.151 and {Hn-i) we see that for 5 small 
enough, \Vk + TVk\aQ < K' , so we can apply the estimates (j5.64p . The new definition of Li{s) 
reflects the third term on the right in the estimates (j5.64p . 

D 

In place of Lemma |5.5| the estimate of substitution errors, we now have: 

Lemma 5.18. 1) For large enough 6q and small enough 5 we have for all k = 0, . . . ,n — 1 and all 
integer s G [0, d] 

(5.71) \elU < C52^^2(.)-i^^^ 

where L2{s) = max(s + ao — 2a + 1, (s — a)+ + 2ao — 2a + 2). 

2) For large enough 6q and small enough 6 we have for all k = 0, . . . , n — 1 and all integer 

s £ [0, d] 

(5.72) {e'Ds+i < C6^ei'^'^-'Ak. 

Proof. Again we use the formulas (I5.25p . By Lemma [53] we have ISg^^Vk + t{I — Sgf,)Vk\ao ^ K' 
for 6 small enough, so we can apply the estimates ()5.64p . When estimating the right sides of ()5.64p 
we use, for example, 

(5.73) \{I - Se,)Vk\s < CSei-'^+K 

D 
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In place of Lemma |5.6| the estimate of accumulated errors, we now have: 
Lemma 5.19. There exist 9q sufficiently large and 5o sufficiently small so that for < 6 < 5o 
(5.74) \En\a < C6^0^^^^^ and {En)a+i < C5^e^^^^\ 

Proof. Since a = 2a — oq, we have ^2(0) > 0, so the proof is the same as that of Lemma 15.61 D 

The new version of Lemma 15.71 the estimate of /„ and g(„, is: 

Lemma 5.20. There exist 9q sufficiently large and 5q sufficiently small so that for s £ [0,q; + 1], 
< 5 < 6q we have 

(a) |/„|. < C520^^W-iA„ 

(b) {gn)s+i < C520^2(s)-i^^ + C0r°"'(5)a+2A„. 

Proof. Since s — a + L2{a) < L2{s), the proof of Lemma [5771 can be repeated here. D 

Induction step. For 6 > sufficiently small, the estimate for the linearized system (j5.63p 
applies to (15. 7j) and gives for s G [0, d]: 

(5.76) |K|s + (K)s+1 < C [\fn\s+l + {gn)s+2 + (l/nUo+l + {9n)ao+2) {\Se^Vn\s+l + {Se^Vn)s+2)] 

Indeed, (j5.67p implies that iov 6 > small enough, Sq^Vu satisfies the requirement (|5.62p r^l For 
the terms involving /„ and gn, except {gn)ao+2, we substitute directly into (15.76P the corresponding 
estimates from Lemma [5.201 For {gn)ao+2 we have 

(5.77) (5„>ao+2 < C (520^2(-o+i)-iA„ + e-''-\g)ao+a+3^ 



where the last term arises from an estimate like (j5.35p with s = qo + 2 and a replaced by a + ao + 3. 
We also use 

(5.78) {Se,yn)s+2 < C50i^+^-")++^ 

and a similar estimate for l^e^l^ls+i, which follow directly from (|5.67p . 
Making these substitutions in (|5.76p gives for s G [0,(5]: 

(5.79) 

\Vn\s + {Vn)s+l < 



c 



For s e [0, a] the parameters oq and a (recall (15.60p ) satisfy: 

(a) L2(s + 1) <s-a 
(5.80) (6) L2(ao + I) + {s + 1 - a)+ + 1 < s - a 

(c) (s + 1 — q)+ < s. 

Thus, we have proved (H„), which is the content of the following lemma. 



Here we use a\ < a. Also, the term \eddiSe^Vn)\ai-i is estimated using the equation (|5.7|) 
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Lemma 5.21 (H„). If 6 > and {g)a+2/^ <ire sufficiently sTuall and 6q sufficiently large, we have 

(5.81) IKU + (V;),+i < (J^-^-^A^ for all integer s G [0, a]. 

Still assuming (H„_i) we now show: 

Lemma 5.22. Suppose n > 1. If 6 > is sufficiently small and 9q sufficiently large, we have 

{a)\C{Vn)\s< 50 n'""'^ for all integer s G [0,d] 

(6) {B{Vn) - g)s+i < 69n'''~^ for all integer s G [0, a]. 

Proof. From ()5.13p we have 

(a) \C{Vn)\s <\{I- Se„_Ji?„-l|. + \en-l\s 

(b) {B{Vn) - g)s+i < ((5,„_, - I)g)s+i + ((/ - Se^_,)En-i)s+i + (e„-i>s+i. 
Using (|5.55p and the above estimates of En-i and e„_i, we find 

\en-l\s<C6^9t^'^''An, 

which imply (j5.82p (a) since oq — a < and L2{s) < s — a. 

The last two terms on the right in (]5.83p (b) are estimated similarly. To finish we use 

(5.85) ((5e„_, - I)g)s+i < CC'"''(5)a ioi s < a - 1, 

and observe that s — a + l<s — a — 1. D 

We now fix 5 and 6*0 as above and check (Hq). 
Lemma 5.23. // {g)a+2 is small enough, then (Hq) holds. 
Proof. Applying the estimate for the linearized system to 

(5.86) ^ ^ . 

B'iO)Vo = Se,g 

we obtain for integer s G [0, a]: 

(5.87) \VoU + (Fo)3+i < C{Se,g)s+2 < C 



9^-(5)a+2, s>a 



^{9)a+2, s <a 

Thus, (Ho) holds if {g)a+2 is small enough. D 

We can now complete the proof of Theorem 15.131 
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Proof of Theorem \5.13[ We have 

n— 1 n— 1 

(5.88) Vn = Vn-l + Vn-l = Vo + ^Vk = ^Vk. 

k=0 k=0 

Let fj, = a — 1. Since 6^ ~ Vk we have by (H„) 

oo oo 

(5.89) 5^|T4U + X^(T4)^+i<<5j^e^2^fc<C7j^A:-i < oo. 

fc=0 k=0 k k 

Thus, for some V as described in Proposition 15. 13( Vfc — t- ^ in E^ and Vk\xa=o -^ 1^1x^=0 in H^~^^ 
(in fact, uniformly for < e < Eq). Lemma |5.2'2I apphed with s = /i — 1 now imphes that ^ is a 
solution of the semilinear system (|5.58p . D 
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A A calculus of singular pseudodifFerential operators 

Here we summarize the parts of the singular calculus constructed in [CGW12] that are needed in 
this article, and we also prove some additional smoothing properties for some of the operators that 
appear in the proof of Propositions 12.41 and 12.61 



A.l Symbols 

Our singular symbols are built from the following sets of classical symbols. 



Definition A.l. Let O C M be an open subset that contains the origin. For m G M we let S^{0) 
denote the class of all functions a : O xW^ x [1, oo) — )• C^, M > 1, such that a is C°° on O xR^ 
and for all compact sets K C O: 

sup sup sup(7' + |e?)-('"-l'^l)/'|9°9g",c7(t;,e',7)l<Ca,.,i^. 
veK g'eRrf 7>i 

Let C^{M.'^ X T), A; G N, denote the space of continuous and bounded functions on M*^ x M that 
are 27r-periodic in their last argument, and whose derivatives up to order k are continuous and 
bounded. 

Definition A. 2 (Singular symbols). Let m G M, n G N, and fix (3 £R'^\0. We let S^ denote the 
family of functions (oe,7)eg(o,i],7>i ^^^^ ^'^c constructed as follows: 

(A.l) y{x',eo,^',k) gM"^ xTxM'^xZ, ae,^{x' ,eo,^',k) = a (eV{x' ,60),^' + —,-f] , 

where a G S^{0) and V G C^(M'^ x T). Below and in the main text we often set 

e 

A. 2 Singular pseudodifferential operators 

To each symbol Os^-y as in (|A.1|) . we associate a singular pseudodifferential operator Op^'^{a) 



whose action on Schwartz class functions u G 5(M'^ x T : C^^) is defined by 

(A.2) Op^>^a)uix',eo):=-^Y. [ e'^^'^'^^''^^ a (eV{x' ,60),^' + —,l) He,k)de, 

where u{S,',k) denotes the Fourier transform at ^' of the A;-th Fourier coefficient of u with respect 
to ^0- When a^^-y is defined as in (jA.ip . below and in the main text of the article, we will often 
write a{eV{x,6Q),X,j) in place of ai;^j{x' ,9q,^' ,k), and an in place of Op'^''^{a). In particular, we 
let Ajy denote the singular Fourier multiplier associated to the function 

A(X,7):=(7^ + |Xp)V2. 

When V{x' ,Xd,6o) depends also on a normal variable Xd > 0, we define the associated family 
of operators depending on the parameter x^ in the obvious way. The pseudodifferential calculus 
takes place only in the tangential directions {x',9q). To discuss mapping properties we first define 
"singular" Sobolev spaces as follows. 
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Definition A. 3. We let 

H''%R'^ xT):= lu 

This space is equipped with the family of norma'^n 



G S'iR'^ X T) : V / (1 + \X\^y |n(e', k)\^d^' < oo I 






|2 ._ 

(27r)« 



Observe that for s fixed the space ff*'^ depends on e with no obvious inclusion if ei < £2- 
However, for fixed e > the norms | • |//s,e ,y^ and | • |//s,e_^2 ^'^s equivalent. 

The next Proposition describes some of the mapping properties of these operators. Detailed 
proofs can be found in ^CGW12] . The constant C is always independent of e € (0, 1] and 7 > 1, 
and we denote the Lp'iW^ x T) norm by | • |o (which corresponds to s = in Definition I A.3p . 

Proposition A. 4 (Mapping properties), (a) Suppose a{eV{x,6Q),X,^) G 5™, where n> d + 1 
and m < 0. Then an ■ L'^{R'^ x T) ^ L'^{W^ x T) and 

(h) Suppose cr(ey(x,6lo),X,7) G 5;j^, where n>d+l andm>0. Then ao : ^'"■^(M'^ x T) ^ 
L2(m«( X T) and 

(c) [Smoothing property] Suppose a{eV{x,6Q),X,'^) G S^"^ , where n > d + 2. Then ud '■ 
L2(R'^ X T) ^ H^^'iR'^ X T) and 

\aDu\H^,e,y < C|n|o. 

(d) Suppose a{sV{x,eo),X,-f) G S^, where n>d + 2. Then an ■ H^^'{R'^ x T) ^ H^^'{R'^ x T) 
and 

Residual operators. We sometimes denote by tq^d an operator that maps L^(]R x T) — ?■ 

L'^[R'^ X T) and satisfies a uniform operator bound 

\ro,Du\o < C\u\o, 

even when tq^d is not necessarily defined by a symbol in some class S^. Similarly, we sometimes 
let r_i^£) denote an operator not necessarily associated to a symbol in S^^ such that 

(A.3) Ir-i^Dulnh^^-y <C\u\o. 

For example, the composition U-i^d ^0,0 = '^-i.D of an operator of order —1 (case (c) in Proposition 
IA.4P with an operator of order (case (a) when tti = 0) is of this latter type. 



^^In this appendix we use j • j instead of (■) in the notation for norms on R** x T, but otherwise we retain notation 
from the main text. 
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Remark A. 5. Observe that a composition of the form rQ^^r^i^D is not necessarily an operator of 
type r_i_£), a fact that is a source of difficulty in the proof of the main linear estimate, Proposition 
12.21 This is the case, for example, if rQjj is the operator of multiplication by V{x' , Oq) G C^(M'^ x T). 
On the other hand we have 

ey(a;',6'o)r_i^Z) = r_i_£, , 

and more generally. Proposition IA.4r d) implies that if o" G 5^, n > d + 2, we have 
(A.4) (Tor-i^D = r-i^D- 

A. 3 Adjoints and products 

In spite of the fact that singular symbols and their derivatives fail to decay in the classical way 
in (^', k, 7), it is possible to construct a crude calculus of singular pseudodifferential operators with 
useful formulas for adjoints and products which, in particular, permit Carding inequalities to be 
proved. This calculus was used repeatedly in the proof of the main linear estimate. Proposition 



In the next proposition a* denotes the conjugate transpose of the M x M matrix valued symbol 
a, while (an)* denotes the adjoint operator for the L^ scalar product. 

Proposition A. 6 (Adjoints). (a) Let a G S^, where n> 2d + 3. Then {(Jd)* — {<^*)d = ?'-i,D- 
(b) Let a G S^, where n > 3d + A. Then {(Jd)* — {<^*)d = '"o.D- 

Proposition A. 7 (Products), (a) Suppose a and r lie in S^, where n >2d + 3. Then 

(^D TD - ((7 t)d = r-i^D- 

(b) Suppose o G S\, T ^ S^ or a ^ S^, r G S^, where n > 3d + 4. Then 

(7 DTD - {(yT)D = ro,D- 

(c) Suppose a G S^"^ , r G S\, where n>3d + A. Then 
(A.5) aoTD- icrT)D = r-i^D- 

Remark A. 8. Observe that when r = t(X, 7) is independent of eV{x,9o) and thus gives rise to 
a Fourier multiplier, the composition (ToT£, = (itt)/) is exact, a fact that has been used several 
times in the proof of Proposition [ 



The equality ()A.5P does not hold in general when a £ S^ and r G S~^, and this is one of the 
main difficulties in the proof of Proposition (|2.6p . There is however one frequency regime in which 
it becomes possible to prove the analogue of equality (jA.SP for (1, —1) compositions, see Proposition 
|A9] below. 

A.4 Extended calculus 

In the proof of Proposition 12.21 and of intermediate results, see in particular the decomposition 
(12. 6p . we use a slight extension of the singular calculus. For given parameters < 5i < 62 < 1, we 



70 



choose a cutoff x^H' j e > t) sudi that 



■/ fc/3 



(A.6) 



< x' < 1 , 



1 on 



suppx^c|(7' + |C?)^/'<52 



kl3 



k(3 



and define a corresponding Fourier multiplier xd in the extended calculus by the formula (|A.2p 
with x'^(C') e )7) in place of (T(ey,X, 7). Composition laws involving such operators are proved in 
|CGW12] . but here we need only the fact that part (a) of Proposition IA.7I holds when either a or 
r is replaced by an extended cutoff x^- 

The following Proposition is essential for treating some of the remainder terms that arise in the 
proof of Proposition 12.41 Part (h) of Proposition IA.9I is used to handle the remainder in the (1, —1) 
composition in (I2.14p (a). and part (c) is used for that in ()2.30p (a). 



Proposition A. 9. (a) Suppose V S C^{W^ x T). Let fiV) be any smooth function, suppose 
a_i(X, 7) is of order —1, and let x^ be a cutoff in the extended calculus satisfying (lA.6p . Then we 
have 

f{V)a-i^DXD = r-i,D- 

(h) Suppose V £ C^{W^ x T) where n > 2d + 3. Let c{eV) he any smooth function and suppose 
ai(X, 7) G S\ and a_i(Ar, 7) G S~^ . Let also x^ be a cutoff in the extended calculus satisfying 
()A.6p . Suppose also that ai is smooth and homogeneous degree one in {X,^). Then we have 

[qi,d, c{eV)] a_i,D xh = r-i,D- 

(c) Suppose V G C^(]R x T) where n > 3d + 5. Let c{eV) he any smooth function and suppose 
ai(X, 7) G S^. Then we have 

[^x'Mi,D,c{eV)]] = ro^D- 

Proof, (a) There is first an obvious L^ bound: 

7 \f{V) a-i,D Xd^Io < 7 \f{V)\L°° |a_i,D Xdu\o < C \u\o . 

Then we need to estimate the singular derivatives {dx^ + f3j dg^/e), j = 0, . . . , d — 1, of the product 
f{V) a-i^D Xd"^- When the singular derivatives is applied to the term O-i^d Xd"") we use the bound 
|X| |a_i(X, 7)1 < C to obtain a uniform L^ bound, and it thus only remains to look at the term 



/3. 



^ £ 



f{V) 



o-i,dXd'"- 



Since V ^ C^, the only difficult term is dg^ f{V)/e. For this term, we use the fact that on the 
support of x^(C'; ~^)7) we have 

|^,7r'<Ce, 

and therefore |a_i^/) Xd'^Io < C'elulo- 

(h) It follows easily from part (a) that 



(A.7) 



d, 



'00 



5-, +/3i— ,c(et/) 



a-i,DXD =r-i,D, j = 0,...,d-l. 
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We now reduce to this case by writing 

d-l d-l 

ai{X,j) = Y,idx,a)Xj + {d^ai)j = J2boA^^ 7) {iXj ) + 60,7 (X, 7) 7- 
i=o j=o 

With 5o = bo J, for example, we have 



&o,D {d,^+f3j^\,cieU) 



a-i,DXD 



[bo,D,cieU)] ( a. , + /3, ^ 1 a_i,i3 Xd + ^o,D 



d.,+(3j^,c{eU) 



a-i,DXD- 



The first term on the right is clearly of type r_i^/5 as a product of an r_i^£) operator on the left by 
an ro,_D operator on the right. By ()A.7p and ()A.4p the same is true of the second term. We leave 
to the reader the verification that the commutator with the last term bo^^(X, 7) 7 gives also rise to 
an r_i^£) remainder (this is even easier). 
(c) The claim follows from the relation 



[V^', [ai,D,c{eU)]] = [ai,fl, {V^>c{eU))] 



and from Proposition IA.7I (b). 



D 



Remark A. 10. Since \eddU\QO,Mo-i < K hy the assumption in Proposition 12. 4^ the result of 
Proposition IA.9I part (a) yields ddQ-i,D = ^-i,D in (|2.14p (a). Similarly, Proposition IA.9I part (b) 
enables us to obtain (|2.14p (b). 

In the proof of Proposition 12.21 we use the following localized Garding inequality for zero order 
operators. As before we write ( = (^'57)- 

Proposition A. 11 (Garding's inequality). Leta{vX) £ S*'(C') andx{v,C) S ^{O) be such that 

Rea{v,C) > c>0 

on a conic neighborhood o/suppx- Provided the corresponding singular symbols lie in S^, n > 2d+2, 
we have 

Re {aD XDU, XDu) > - \xdu\1 \u\l. 

2 7 
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B An example derived from the Euler equations 

In this appendix we explain in a particular example how one can derive a single nonlocal nonlinear 
equation that governs the evolution of the amplitude function a, which itself determines the leading 
profile V*^, see Proposition 11.231 In the process we provide explicit constructions of a number of 
the objects that appeared in our earlier discussion of approximate solutions. 

As in |CG10) . we consider the linearized Euler equations in two space dimensions to which we 
add a nonlinear zero order term (we slightly change notation compared with the introduction). 
More precisely, we consider the following system 

'dtV' + Ai d,,v' + A2 d,,v' + D(y^ v') = o, (t, xi, xs) g ] - oo, r] x m^ , 



B V'\^,=o + *(y^ V')\^,=o = e^ G{t, xuMt, xi)/e) , {t, xi) G ] - 00, T] x 
t<o 



(B.l) 

where the 3x3 matrices Ai,A2 are given by 






— V 





c^v 


















u 
Ai := I -c'/v I , ^2 := I u 

c^/v 

and the parameters v, u, c are chosen so that 

t;>0, <u < c. 

The latter assumption corresponds to the linearization of the Euler equations at a given specific 
volume V with corresponding sound speed c, and a subsonic incoming velocity (0, u) (observe the 
difference with [CGIO] ). We also assume that D in (jB.ip is a symmetric bilinear operator from 
R^ X R^ into R^, and that ^ is a bilinear operator from R^ x R'^ into R^ (why we choose R^ is 
explained below). 

For such parameters, the operator dt + Ai dx^ + A2 9^2 iii (jB.ip is strictly hyperbolic with three 
characteristic speeds 



Ai(ei,6) := n6 -c^ei+^i , A2(6,6) := ^6 , A3(ei,6) := «6 + c^ef +ei • 

There are two incoming characteristics and one outgoing characteristic, so B should be a 2 x 3 
matrix of maximal rank. The choice of B is made precise below. Of course, the source term G in 
(IB.lj) is valued in R^. We assume moreover that G is 1-periodic and has mean zero with respect to 
its third variable ^o- We choose a planar phase (pQ for the oscillations of the boundary source term 
in (iBljl : 

(po(t,xi) :=Tt + rixi, (z, ;?/) / (0, 0) . 

The hyperbolic region Ti can be explicitly computed and is given by 

7^ = {(r,?7) gRxR/|t| > Vc2-n2|r?|| . 

For concreteness, we fix from now on the parameters (r, 77) such that rj > and T = cr]. In this 
way, we have (r, 77) G Ti. 
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We determine the planar characteristic phases whose trace on {x2 = 0} equals (pQ. This amounts 
to finding the roots uj of the dispersion relation 



det 



Tl + r]Ai + ojA2 



0. 



We obtain three real roots that are given by 

2M 



Wi 



■r], ^2 := , u 



3 •- 



— 17, M :-- 
M- 



G]0,1[. 



l-W^- - - ivi- c 

The associated (real) phases are (jjiit^x) := 4>Q{t,xi) +c<JjX2, i = 1,2,3. The relations 

T + Ai(?7,a;i) = Z + >^i{rii^2) =Z + A2(r/,W3) = 0, 
yield the group velocity Vj associated with each phase (/>j: 

l-M^ f-c\ f-c\ /O' 



vi := 



1 + M2 



V2 := 



u 



V3 



u 



Hence the phase (j)i is outgoing while (j)2,4'3 are incoming. With the notation of the introduction, 
we can also compute 

1-M2 



ri := 




l-M^ 



' 2(1 + M2) 



1/v 
-1/c 

1/u 



r2 ■■ 




rs :-- 



h:-- 



1 



1 + M2 




from which one can obtain the expression of the projectors Pi, P2) -fs as well as the expression of the 
partial inverses Ri, R2, Rs- The stable subspace at the frequency (r , rj) is spanned by the vectors 
r2 , rs . The matrix B in (IB.ip is chosen as 



5 := 



V 
u V 



so that we can choose e := r2 — rs as the vector that spans keri? n W{t,7]). The reader can 
check that all our weak stability assumptions are satisfied with this particular choice of boundary 
conditions. (We skip the details that are just slightly more complicated than those in [CGIO] .) The 
one-dimensional space -BE'^(t, r/) can be written as the orthogonal of the vector b := (n, — c) . 
The leading profile V'' and the corrector V^ satisfy 

V° = VL = a2{t, X, 02) r2 + a^it, x, 0s) rg , V^^ = n{t, x, 0i) n . 

Moreover, there holds 

V°(t,xi, 0,6*0, 6*0, 6*0) = a{t,xi,9o)e = a{t, xi,0o) {r2 - rs) , 
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where the scalar function a is 1-periodic with respect to ^o and has mean 0. The Fourier coefficients 
of a are denoted ak, k £ Z, where oq equals for all time t. Since the functions o"2, cr^ satisfy the 



transport equationO 

dt(72 + V2 • Va;C72 = dta^ + V3 • V ^(73 = , 

and vanish for i < 0, we obtain the expressions 

t ,xi + —-,92) , a3{t,x,93) = -a[t ,xi,93 

u M J \ u 

To compute V^^i, we must solve 

(B.3) Eout (Lid) VL + ^2"' D(VL, VD) = , (here E^ut = E,) , 

and we thus need to determine the resonances between the phases. A simple calculation shows that 
there is a nontrivial n G Z^ satisfying ni (pi = n2 (/>2 + ^3 <p3 if and only if M^ is a rational number. 
We thus assume this to be the case from now on. The resonance between the phases reads 

ni:=q, n2 := p + q , ns := -p , with ^ _ j^2 = ~ ^ 

and it is understood that p, q are both positive and have no common divisor (for instance p = q = 1 
when M equals l/\/3). Expanding the quadratic term ^{Vinc^^inc) ™ Fourier series, and using 
the relation 

Ci = Z U Z n2 
we obtain (use the expressions ()B.2p ) 

El {A2' D(VL, VD) = -2 E «^(P+.) (* - ?' -1 + S) --^P (* - ?' -0 «'^^'''^ 

kez 

Pi^^iD(r2,r3). 
In terms of the interaction integral, we obtain the expression 

El (^2iD(vL,vL)) 

= -2 f\a)^ft-^,xi + ^,^ei-^e3) a(t-^,xue3) d93PiA^'B{r2,r3), 
Jo \ u M n2 n2 J \ u / 

where (a)„2 still denotes the action of a under the preparation map that retains only Fourier 
coefficients that are multiples of 712 . Consequently ()B.3P reads 

. , /^ 1-M2 ^ 1-M2 ^ \ 

(B.4) Ot K COx, ^ UOxo Ti 

^ ' V 1 + M2 "^1 1 + M2 ^V 



d / (a)„2 [t 



X2 X2 ni '^3^] f. X2 . 

— '^1 + ITF' — ^1 ^3 a[t ,xi,93] d93, 

u M n2 712 / V n 



*^ Observe that there is no zero order term in the transport equations because the zero order term in ()B.1[I has only 
a quadratic part. This choice has been made for the sake of simplicity. 
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with 

X — m"^ 

The transport equation ()B.4p is solved by integrating along the characteristics and we obtain the 
expression 

The Fourier series expansion of ti reads 

r. . , , V- /"* f2s-(l-M^)t ^ l-M^ , A 

Ti(i,xi,0,ei) = d^ / afc(p+,)( ——2 ,xi + 2c——^{t-s)\ 

/2s-(l-M2)t 1-M^, A , ,... 



The equation governing the amplitude a reads 

b-{{a^r^{e,e)+n\^,=oBn-BR{L{d)Vl,)\^,=o)=b-G, 

where functions are evaluated at X2 = and Oi = 02 = 9^ = 9q. Since we already have the 
expression of ri in terms of a, the only task left is to compute the trace of the term B R {L{d) Vf„^). 
Recalling that i?2 ^2 = ^3 '"s = 0, we have 

BRiL{d)Vl,)U=o = {BR2A^\2 + BR3A^\3)dta+iBR2A^^Air2 + BR3A^^Air3)d,,a, 

with a the unique primitive function of a with zero mean. Using the expressions of R2, R3 in terms 
of the projectors Pi, P2, P3, which themselves can be obtained from the vectors ri,ii, we get 



b.iBR2 A-,^ r2 + BR3 A~,' n) = -^^^^%^ 



h.{BR2A-,'A,r2 + BR3A,'A,r3) = ^^^^^^^^^. 

The fact that both quantities are proportional one to the other with a factor — c comes from a 
general fact, see [CGlOl Lemma 5.1]. 

The function a should therefore satisfy the amplitude equation 

y^ '- [dta - cd,,a) + b ■ ^(e, e) (a^)* + b-Bn tiU.,=o = b-G, 

or equivalently 

(B.6) "j^ ^ {dta - cd,,a) + b ■ *(e, e) dg, (a^) + b-Br, dg,T,\,^=o = b ■ de,G . 
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Let us define the two constants 

"^ ■= uv{i+~M^) ^ ■ ''(^' '^ ' "^ ■= TTJ^ ^' • ^^ °("^' ^^^ • 

Then equation (1B.6P reads 

9ia - cSxiO + ai de^ia ) + a^ S^o -71x2=0 = ,. ~..o^ b ■ de^G , 

a uv [1 + iW^j 

where the derivative dggTi/d\x2=o is computed from the relation (|B.5p : 

. n ni f' f\ ( 2s-{l-M^)t 1 - M^ m ns 
deo-7\x2=a = — / / Oe^a „2 2 '^1 + ^*^ 1 I /i#2 (^-g . — ^0 

d "2 Jo JO V 1 + ^ 1 + ^ "2 "2 

In terms of the Fourier coefficients at , the latter equation is seen to be equivalent to the infinite 
system of transport equations 

dtak-cdx^ak + 2i'Kkai > a^' a^-fc' = 2 z vr A; ~n^9\ b-Gk, k^qZ, 

k'&L ^ ' 

and 

dtakq - cdxj^akg + 2i-Kkqai ^ Ofc' akq-k' 

fc'ez 
o. , /"* /2s-(l-M^)t l-M^ , A 



We recall that the coefficient oq vanishes. 

In the special case M = l/\/3; the above system reduces to 



dtak - cdxiak + 2i-Kkai ^ Ofc' afc_fc' 

fc'ez 

/■* fis-t , \ (Zs-t c , A , 

+ 2i7:ka2 a2k { — 5 — ,xi + c[t- s) j a-t ( — 5 — ,3:^1 + 2 (* ~ ^) ) '^^ 



rj 
2iiT k —=— b ■ Gk , k £ Z , 
Auv 



with parameters ai,a2 computed from the nonlinearities D,^ in (IB.ip : 

ai := -^b-^(e,e) , a2 := ucrjii ■ Ar, D(r2,r3). 
4:UV — 
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